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Am 1. 6. 1960 erlag o. Prof. Dipl.-Ing. Dr. techn. Dr. Ing. E. h. Ernst Chwalla im 
Alter von 59 Jahren véllig unerwartet einem Herzschlag. 


Prof. Chwalla entstammte einer alten Wiener Familie. Schon als Student des Bau- 
ingenieurwesens an der Technischen Hochschule Wien fiel er durch seine hervorragen- 
den Leistungen auf. Er legte simtliche Priifungen von der Matura bis zum Rigorosum 
mit Auszeichnung ab und war dann als Assistent an der Lehrkanzel fiir Briickenbau 


Ree 


bei Prof. Hartmann tatig. Hier habilitierte er sich im Jahre 1928 fiir das Fachgebiet 
Baustatik mit der Habilitationsschrift ,,Beitrage zur Stabilititstheorie“’. Diese Schrift 
kennzeichnete bereits das kiinftige Arbeitsgebiet Chwallas, auf dem er so Hervorragen- 
des leisten sollte. Es folgte nun eine rege Publikationstatigkeit, und schon im Alter von 
28 Jahren erreichte ihn die primo et unico loco-Berufung als Vorstand des Instituts 
fiir Baumechanik an der Deutschen Technischen Hochschule in Briinn. Er wurde wegen 
seines jugendlichen Alters zunachst als Extraordinarius bestellt. Im Jahre 1937 erfolgte 
dann seine Ernennung zum Ordinarius. 

Die Zeit bis 1945 war ausgefiillt mit einer intensiven Lehr- und Forschungstatigkeit, 
als deren Frucht nicht weniger als 61 wissenschaftliche Arbeiten erschienen. 1943 wurde 
Chwalla zum ordentlichen Mitglied der Deutschen Akademie der Wissenschaften in 
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Prag gewahlt. Schon seit 1939 wirkte er mafgeblich an der Schaffung der deutschen 
und spater auch der ésterreichischen Knicknorm mit; die heutige DIN 4114 und die 
O-Norm 4300/4 sowie ihr in Kiirze erscheinender Neuentwurf sind in wesentlichen 
Teilen sein Werk. Daneben war Chwalla als Experte der Industrie und der Behorden 
an den Entwurfsarbeiten und Festigkeitsuntersuchungen fiir verschiedene grofe Bau- 
vorhaben in Deutschland beteiligt. Erwahnt seien hier nur die Elbe-Hochbriicke in 
Hamburg mit 750m weit gespanntem Mittelfeld und die 240 m hohe Rahmenkuppel 
des Miinchner Hauptbahnhofes. 

Der Zusammenbruch Deutschlands im Frithjahr 1945 brachte auch das Ende der 
Deutschen Technischen Hochschule Briinn. Chwalla verlor nicht nur seine Stellung 
als Ordinarius, sondern auch sein gesamtes in Briinn befindliches Vermégen, alle Manu- 
skripte und seine groBe Privatbibliothek. Mit ungebrochener Energie ging er jedoch 
daran, nun in Osterreich sein Leben neu aufzubauen. Er entfaltete eine auBerordentlich 
weitgespannte Tatigkeit als Ingenieur-K onsulent fiir die Vereinigte Osterreichische Eisen- 
und Stahlwerke-A. G., fiir die Vorarlberger Illwerke A. G., die Oberdsterreichische 
Kraftwerke A. G., die Niederédsterreichische Elektrizitatswerke A. G. und die Siemens- 
Schuckert-Werke. Die Tatigkeit Chwallas in dieser Zeit ist mit der Entwicklung der 
dsterreichischen Wasserkraftwerke auf das engste verkniipft. Er war Amtssachver- 
stindiger der Obersten Wasserrechtsbehérde beim Entwurf, der Berechnung und der 
Uberpriifung von Talsperren und damit mageblich verantwortlich fiir den statischen 
Entwurf fast jeder der groBen Staumauern Osterreichs. Nebenher liefen noch seine 
Arbeiten in den verschiedenen Fach- und Normenausschiissen des In- und Auslandes. 


Unmittelbar nach seiner Riickkehr nach Osterreich im Jahre 1945 war Chwalla 
zum Honorardozenten fiir Stahlbetonbau an der Hochschule fiir Bodenkultur ernannt 
worden. In der Folge versuchten zahlreiche auslandische Hochschulen, ihn durch 
primo loco-Berufungen fiir sich zu gewinnen. Chwalla folgte diesen finanziell oft auBer- 
ordentlich verlockenden Angeboten nicht, sondern nahm im Jahre 1955 eine Berufung 
als Ordinarius und Vorstand des Instituts fiir Baustatik an der Technischen Hochschule 
Graz an. Im Jahre 1954 erhielt er das Ehrendoktorat der Technischen Universitat 
Berlin-Charlottenburg, und im Jahre 1959 wurde er zum ordentlichen Mitglied der 
Osterreichischen Akademie der Wissenschaften gewahlt. 


Chwalla hat gegen 100 wissenschaftliche Arbeiten aus der Baustatik, der Elastizi- 
tats- und der Plastizitatstheorie verdffentlicht. Er besaB in ungewéhnlichem Mae die 
Fahigkeit, das Wesentliche eines Problems sofort zu erkennen, er besa aber auch die 
unermidliche Zahigkeit, die Lésung bis zum numerischen Ergebnis voranzutreiben. 
Fast alle seine Arbeiten bedeuten Vorst6Be in unerforschtes Neuland. Vor allem die 
Stabilitatstheorie verdankt ihm auRerordentliche Fortschritte. Sein Vortrag war von 
einer unerhorten Klarheit, witzig und mitreiBend. 

Chwallas friiher Tod trifft seine Freunde hart; fiir die 6sterreichische Bauwissen- 


schaft aber bedeutet er einen Verlust, der noch auf viele Jahre hinaus fiihlbar bleiben 
wird. H. Parkus 


Compressibility Effects in Circumferential Inlet Distortion in Axial 
Compressors* 


By M. Z. vy. Krzywobloecki, Urbana, Illinois 
Part bt 
With 1 Figure 


Summary. This is the second paper of a series of papers dealing with the compressibility effects 
in circumferential inlet distortion in axial compressors. The first paper referred tJ the solution 
of Ehrich valid in the incompressible fluid domain. In the present paper a technique of extending 
the proposition of Yeh to the compressible fluid flow regime is demonstrated. The two examples 
worked out show quite readily that the generalization is a straightforward one and the amount 
of the additional calculations due to the compressibility effects is small. 


Introduction 


The present paper is the second one in the series of papers dealing with the 
compressibility effects in the circumferential inlet distortion in axial compressors. 
In the first paper of the series the author generalized to the compressible flow 
regime the method of Ehrich of analysing the circumferential inlet distortion 
phenomena in axial compressors in the incompressible fluid flow domains. In it 
the author derived some similarity rules which enable one to generalize the re- 
sults from the incompressible flow domain to compressible flow regime. Besides 
that a straightforward generalization to the compressible flow was worked out. 
This technique is applied in the present paper to the approach to the circum- 
ferential inlet distortion proposed by Yeht. 


1. Yeh’s Solution in the Incompressible Fluid Flow 
1.1. Basic Equations 


Yeh considers the case of a plane cascade of infinite number of blades and of a finite 
span between two walls. The cascade is approximated by the actuator disc located at x = 0. 
The coordinate axes are located as follows: x—axial direction, y— circumferential 
and z—spanwise direction. The undisturbed velocity vector has the components 
(U, V, 0), the absolute main flow angle is «=tn-!(V U~), the velocity of the rotor 
blade is equal to k U, with k=q-', @ (flow coefficient) = U (w n)~, the relative main 
flow angle is B=tn- [(V +4 U)U~ ]. The main flow is uniform along the span, 


'* The results presented in the present paper were obtained during the author’s association 
with the Compressor Aerodynamic Design and Development Unit, General Electric Company, 
Cincinnati, Ohio, U.S.A., with Mr. R. A. Novak as its supervisor. As in the past, the author wishes 
to express his deep appreciation to Mr. R. A. Novak for the numerous discussions on the subject 
in question. 

** Part I: Osterr. Ing.-Archiv 13, 214 (1959). 

1 Yeh, H.: An Actuator Dise Analysis of Inlet Distortions and Rotating Stall in Axial Flow 
Turbomachines. General Electric Co., Cincinnati, Ohio, Aircraft Gas Turbine Division, Technical 
Information Series, No. R 57 AGT 628, Sept. 1957. 
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the blades are untwisted, the region upstream of the cascade is denoted by the sub- 
script 1, and that downstream of the cascade by the subscript 2. Yeh represents the 
inlet distortion by a velocity perturbation consisting of two parts: (i) the distortion 
velocity components w and v, carried along by the main velocity: Wp (x, y, 2)= 
= uy ( tna— y, 2); vo (x, ¥, 2) =U (w tna—y, z); and (ii) the induced velocity com- 


ponents wu’ (x, y, 2), v' (a, y, z), w’ (x, y, 2). Yeh finds the induced velocity-field, which. 


yields the exit distortion uw), and vy), for the given inlet distortion uw, and Voi LBC 
used equations are the linearized Euler equations and the continuity equation for an 
incompressible flow: 


UU! + VU, y=— OP sx; (ere 
Ov',.+ Vv',y=— 07 Psy; (12?) 
Uw',.+ Vw',y=— 0 pez; (iL3) 

(1.1.4) 


Wee yw, 7=0. 


The inlet distortion is represented by a velocity perturbation u,, along x- and vp, 
along y-axis far upstream at x= — oo. Due to this inlet distortion, the velocity far 
downstream at x=-+ oo is perturbed by the amounts w,, and vp». 

Thus one has: 

Far upstream: U +91 (— ~, ¥, 2); Vi tvo1 (— ©, ¥, 2); 95 (1.1.5a) 


Far downstream : U + U2 (+ -, y, 2); VetVo2(+ 0, y, 2); 9. (1.1.5b) 


When p’ is cancelled among any pair of the first 3 equations, one obtains 3 equations 
of the form: 


U (Uy — 0 x)ox bY (tt’,y— 0 ox)59 = 0; ete., (1.1.6) 


which means that the vorticity components are equal to zero or that there exists the 
velocity potential ®’, i. e., ®’,,=w’, etc. Then the continuity equation yields a Laplace 
equation for ®’ whose solution is proposed by Yeh in the form: 

Upstream of the cascade: 


© = 2 D) Anum XP [(An? + uz) x] exp (iAn Y) COS Mm 23 (ga 
Downstream: wes 
®,/ = D>) DS) Bimexp [— (An? + Mm) @] exp (tA, y) COS Lm 2; (1.1.8) 
where: ae 
Ai OT in = 90 I Leet, (es 0 Le Dee (121.9) 


with L, denoting the wave lenght along y, L, denoting the blade span. The coefficients 
Aym and Bym are, in general, complex numbers. 
Yeh finds: 


een Ale (An? + wi)” exp [(An?-+ wz) 2] + exp (% An Y) COS Mm 25 (1.1.10) 


v1! = D>) dD) Anmt dn XP [(An® + uz)” 2] - exp (t An Y) COS bm 23 GS eg) 
Di rer ae eas ais exp [(A,?-+ u2)” x]: exp (i An y) SIN Um Z, (1.1.22) 


with similar equations for w,’, v9’, we’. 


As the matching conditions Yeh assumes the following conditions: 
(i) the continuity equation at the disc; 
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(ii) the condition that the spanwise velocity should undergo no change across the 
actuator disc; 

(ili) the magnitude of the exit flow angle £,; 

(iv) the loss of the total pressure across the cascade. 


Below, we shall consider the modification of this proposition to compressible fluids. 


2. Compressible Flow 
2.1. Fundamental Equations and Similarity Coefficients 
Assume that the density in the axial compressor flow, denoted by the symbol 9, 
consists of three parts: o,= e+ 0)+ 0’, where the quantity o refers to the incoming 
undisturbed flow, 0 is due to the flow distortion and 0’ refers to the induced disturbances. 
The pressure, p,, is assumed to be of the similar structure, i. e., p= p+ p +p’, with 
p being a constant, independent of the coordinates {x, y, z}. The i incening flow (o, U, V) 
is steady, inion and is not a function of the coordinate system {x, y, z}. Then the 
Kuler equations of motion are: 


eo. (U+u,+ wu") (OU +u,tw’).+(Vtuy+v’) (U+u,4+u’'),,  w(U +u+wv),J= 
=—(P+ PtP). ter Ps (2.1.1) 

oa ((U +u,+ 4’) (V+ 4-4 ed L(V+e%t2 Uw) A Ve 0), 9+ wv’ (V+ 05+’), .]= 
ae Pete aOR oo (2.1.2) 


Or (Oty tu’) wet (V+m+0') wy tw w',.J=—(PtPot+P')2c+aP:, (2.1.3) 
with {P,, P,, P:} denoting the external force components; the continuity equation: 


(0: U1), + (o Vi)sy + (Oe W,),.=9, (2.1.4) 
where: 


U,=U+utw; Vi=Vtontv; Wi=w'. (2.1.4a) 


Introduce a linearization into the system of equations with an approximate assumption 
that the pressure p’ and force contribute to the variation in the {u’, v', w'} — com- 
ponents only. With u,+ wu’ <U; vy9+v’ <V, or even u+u' <U, ete., and a,+ 0’ <e 
the equations of motion take the form: 


OC Bo Vie) =p OP je Ewe Zs ( 
o (Uw'",.+ Vo») = Beery ene (2.1.6) 
o(Uw'",.+ Vu',,)=—p ztef,., ( 


F (x, y, z) being the force potential. Differentiating equations (2.1.5) and (2.1.6) with 
respect to y and 2, respectively, (2.1.6) and (2.1.7) with respect to z and y, respectively, 
furnishes three equations of the form: 


o [U (w',y—v',x)oet V (w’,y—'sx),y] =9, ete., (2.1.8) 
which shows the existence of a potential function ®’ (x, y, z), ®’,,=w’', etc. The equa- 
tions of motion can be written in the forms: 


LU ; dV 
a = — (PotD) «HOF ss OG =— Pot P's + OF», ete. ; (2.1.9) 


let us multiply the left hand sides of these equations by U, dt, V,dt, etc., respectively, 
and the right hand sides by U, dt = (da/dt) dt dx, dy, etc., respectively. The addi- 
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tion of these equations, so transformed, furnishes the Newton equation taken along 
a streamline: 


a. (U,dU,+ VidVit W,dW,)=ao.d [4 (U? t+ VP W2)\= 


=ad(sP)=a9 dg=—dp+odk. (2.1.10) 
From the definition of the velocity of sound, (d p,/d @:)s=c?, one has ¢d 9, =c~* dp, OF 
Cee O1 [ — (U, Ure ViVies —- W, Wx) + Poly etc. (2.1.11) 


6 


Assume a running coordinate, “‘s’’, along a streamline, with the origin (s = 0) located 
in the plane of the actuator disc (— upstream, + downstream); inside a rotor there 
takes place a heat addition to a particle of the medium moving along a streamline; 
moreover, a certain amount of work is done by the external force. Let the symbol 
Q (s) denote the total heat amount added to a particle from s=0 up to the point in 
question, and F' (s) 1h, P (s) ds, the total work of the force P on this distance. Then 


the equation of energy along a streamline has the form: 


b (U2 + V2) + (y—1) 4c =4 (UP + VAt We) + (y—1) 402 — 


— [Q (s) +F (s)]; (2.1.12) 
with (U?+ V?— U2— V?,— W,?) = 0, one may assume: 
C2 w ¢2= 2+ (y—1) [Q (s) + F (s)]. (2.1.13) 


In this approach Q (s) and F (s) may be taken to be some average constant magnitudes, 
or parametric functions not being subject to differentiation or integration operations. 
Inserting eq. (2.1.11) into the continuity equation (2.1.4) joimtly with eq. (2.1.13), 
and applying the linearization process analogous to the classical one, furnishes the 
formula: 


[1 — (U + eg)? 637] ®,2e + [L—(V + 09)? 637] P's + O22 + 
+0; [(U +9) Pe + (V +0) Py] =0. (2.1.14) 


Unless the force components are very large, one can neglect the last term. The following 
notation will be used: 


[1—(U + u,)? c7?7]= 6.7; [1—(V + v9)? 72] =8.?. (2.1.15) 
Let us introduce the transformation of coordinates {x, y, z} > {Xin, Yin, Zin} and 
D (x,y, 2) > ©5, (Wns Yin Zin) with (XH, Yn, Zn) and Ot (%%; Yin, 2m) Cenoting the 


coordinates and the velocity potential function in the incompressible fluid flow domain. 
Let: 


Aggie Os hy a US Aen = 8s Ae Pal Ti, Gate =O (aaa) (2.1.16) 


where /,; 4,; 4,; Aw’ are some constants to be determined. The transformed eq. (2.1.14) 
takes the form: 


Bae ho Hk Dj, 2 Xin Xin a Bo" he’ i PD, 2 Vin din = ho Me Di, » Zin Zin = 0 S (2. 1 ED 
We require that all the coefficients in eq. (2.1.17) be equal one to each other: 
Bu® ALA = Bo A? = 22". (2.1.18) 


As the next condition we take the boundary conditions at the surface (i. e., «= 4%,, 
Y= ys): 

dys F i , 

de, = (V tote’) (U +9 pes)! (V +09) (U + rug) + 0! (U 49) 1 = 


= (V+%) (T+ tig) A+ (T +040) (yer ge (2.1.19) 
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Assume the relations: 


(U + U9) = Au (UO + Up)ins (V +09) = do (V + 0 9)in- (2.1.20) 


Then one has: 


dy./dx, = (Ay AYins) | (Ax A2ins) ia [A (V aie Vo)in |/ [Au (U + Molen =e 


+ (Ag’/Ay) [Au (U + Ug )in]=? Bin s yin (2.1.21) 
We require again that all the coefficients in eq. (2.1.21) be equal one to each other: 
Ay | Ax = dol hu = Ae! |(Ay du) - (2.1.22) 
Thus we have four equations [(eqs. (2.1.18) and (2.1.22)] in six unknowns: A,; Ay; 2; 
ay eae ae 
Assume two values: 
Aga MA Sas (22.23) 


and calculate the four other values: 
Aes Uy == Bp Ave orale ee bee (2.1.24) 


thus obtaining the following similarity rules: 


tu! =O". = (B2/B.) Pi, sg, = (Be?|Bu) Ui (2.1.25a) 

v' =B',, = (Be?IBo) Bossy = Br U3 (2.1.25) 

wi =O 5 = (8,711,) Ora = Be Wan (2.1.25¢) 

D(x, y, 2) =P," ©;, (Lin, Yin, Zn) - (2.1.25) 
Similarly one gets: 

(U + U9) = Bu (U + Up)in3 (V + 9) = Bo (V + U)in; (2.1.26a) 

C= ae = pe ans be = (2.1.26b) 


The value of the density, o,, is calculated from the equation of energy with the 
use of: 


p: p= (@: 0-1) exp [(S,— 8) c71], S=entropy, (2.1.27) 
=O (y—l yh (U8 V2 (OC bea’)? = (V + 09+ 0')2 = w'?) 
+@ (s)+F (s)] @ p-?+1} exp [—(S,—S) eF1]}/0-, (2.1.28) 


where the functions Q (s), F (s), S;, S are known and given. The pressure can be cal- 
culated from eq. (2.1.27) and the temperature from equation of state. 


2.2. Straight Forward Approach 


Assume the existence of the velocity potential ’ =@’ (x, y, z) (see the derivation 
below eq. 2.1.8) which possesses the form chosen by Yeh, i. e., it satisfies the conti- 
nuity equations of an incompressible fluid flow and the Laplace equation which ori- 
ginates from the continuity equation of an incompressible fluid: this is for — co <a <0: 


@,' = YD Anmexp [(An? + m2)” x] exp (¢ An y) COS (Um 2); (2.2: 
Oe <= Cols 


@,'= > S| Bam Xp [— (An? + un?)” az] exp (2 An y) COS (fim 2); (2.2.2) 


n m 
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where: 
An=2n2 L51; L,=wave lenght along y,  (#=20, 1. 2, ee ee) 
Um =2mnL;1; L,=blade span, (m=0, 1,2, ...). (2.2.4) 


The coefficients A,,, and B,,, are in general complex numbers. Consider the conti- 
nuity equation of a compressible fluid flow in the steady state conditions: 


(01,x OF Ory Vi+e.,< W:)+ (Cie y gar Wr2)=0. (2.2.5) 


Let us choose the values of the velocity components in the forms identical to those 
in the incompressible fluid flow domain. 
From (1.1.7) and (1.1.8) we have for — 0 <a <0: 


Uy == pelea a*nm (An? + bm")? exp [(An” “pn x | exp (2 An y) cos (Min z); (2.2.6 a) 


vy! = YY Aum (6 An) EXP [a2 + 2) 2] exp (i Aan 9) 008 (ptm 2); (2.2.6) 
wy! = — YY Am fm OXP [(An? + fm?) ? %] exp (0 An Y) SiN (Um 2); (2.2.6c) 


and for 0 <#% <-+ oo: 


tbe’ = D>) DS) Bam (— (An? tpn?) *) exp [— (An? +fn®)® #] exp(t- An Y) COS. fm 25 (2.2.74) 


Ve) = 2) 3) Bum (t An) EXP [— (An? + Um?) %] exp (6 An Y) COB pm 25 (2.2.7b) 

wo’ = — 3) >) Bam bm XP | — (An? + pm?) &] exp (1 An Y) SIN Mn 2; (2.2.76) 
or ie 

Vie 0 aie Abi os, “Uge= ==) Varn thas (2.2.8) 

Deyo == hey (An? + fm?) ~” ; eae ne (2.2.9) 

Teme (220 1)? (im L1))-2?; 1, m—=0, 1, 2). 5 te blade span a2 oe 


The velocity components U and V are constant; the velocity components w’, v’, w’, 
satisfy the continuity equation of an incompressible fluid flow, 1.e.: 


Web Uy Wz = 0. (2.2.11) 


Consequently, we assume that the velocity components uw» and v, in the form assumed 
yee lac = 
Ug = Uo (x tan «— Y, Z); Vo= Vp, (x tan « — Y, Z), (2.2.12) 


should also satisfy the continuity equation of an incompressible fluid flow, i.e.: 


Uo,x (wv tan «— y, 2) + 4, (w tan «— y, 2) =0. (2:2.13) 
They may be assumed to be of the value: (— o <x <0): 

Uj, =C exp (+atana+y) f, (2); (2.2.14) 

Vo1 = Uo, tan a= C tan «exp (+atana+y) f, (2), (2:2.15) 


where /, (z) is any arbitrary function of z. 
In an analogous way one can write in [0 <a <+ oo}: 


Ugg= D exp (— x tan a+ y) f, (z); (2.2.14a) 
Vp —D tan « exp (— x tan «+ y) fy (2). (2.2.15 a) 
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Other possible solutions of this equation are: 
exp [i (— x tan a+ y)]; sin (x tan « — y); cos (x tan «— y), 


or any sum of these expressions. 


This implies that the second term in eq. (2.2.5) is identically satisfied by the choice 
of the forms of the velocity components. We assume that the density can be separated 
into three parts: @; 0); 0’. Let us assume now that the part 0, is correlated only with 
the velocity components u, and vy but not with wu’; v’; w’; and that the part 0’ is corre- 
lated only with the velocity components wu’; v’; w’; but not with uo; v». Eq. (2.2.5) 
reduces to the following two forms, with o =const.: 


Qo,x (U ++ u’)+ 00, (V+ v9 +0’) + 00, vw’ =0; (2.2.16) 
o’.x(U+utw)+o',,(V+t%+’)+ 0',,w'=0. BA) 


A solution of a differential equation depends upon the coefficients in the terms of this 
differential equation. This means that a solution 0, of eq. (2.2.16) cannot depend upon 
the (u’, v’, w’)— components and they must be neglected. Analogously, in eq. (2.2.17) 
we have to neglect the terms (U+u,; V+v,). The reduced equation (2.2.16) with 
V=U tan « and v,=u, tan « [see eqs. (2.2.1) and (2.2.15)] has the form: 


00,x (U + Uo) + @0,, tan « (U + u9) = 0. (2.2.18) 
This has a solution of the form: for — «o <2 <0: 
0, =E [exp (x tana+y)]g (z), (2.2.19) 
where g (z) is an arbitrary function of z; for 0 <<a <-+ oo: 
0o=# [exp (— x tan+ y)] 9 (2). (2.2.20) 


Comparing these solutions with solutions for w), v9, eqs. (2.2.14) and (2.2.15) one can 
write: — 0 <4 <0: 

09 = U 0; By U1; (2.2.21) 
0<% <+ co: 

Qo, = U 02 Be Ugo; (2.2.22) 


where the coefficients H,, H, should be determined experimentally. One may also 
propose the solutions in the form: 


Cn =F, 0,0 Wis Cor %1 =F a U; (2.2.13) 
; Ops Ugg 15 U- (2.2.24) 
Thus, having wu; (¢=1, 2) we may calculate 0 ;, and having og); we can calculate 
i 5d Oia Al 1 
Consider the reduced equation (2.2.17): 
On’ +o0',v +0',,wW =0, (2.2.25) 
or in [— ~<a <0]: 
01,2 + O1,y t Vam— 01,2 Kam tan [dm z= 9, (2.2.26) 
and in [0 <%¥ <+ oo]: 


02,2 — 02,y © ¥nm + 02,2 Knm tan Um z=0. (222,27) 
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Solutions are of the forms: 
[—-coo<a <0]: 


o,) = D>) Dd) Bu lexp (2 e—1 (v5; y)]} (sin ez) Vem; 2.2.28) 
5 t 

Toe ae {exp [—2 a4 -+2 (07) y)]} (sin fw z)* (Kuwteo) | (2;2::20) 

Consider the neglected terms in eq. (2.2.16): 
00,x & + 00, V' + 00,, wv’ =0. (2.2.30) 
This equation is mei to eq. (2.2.25) and consequently its solution is of the form 
similar to eq. (2.2.28), i.e., it depends upon the wave lengths L, and L,. Since we 
have assumed above that aie part 0, depends only upon the velocity components wu, 
and v, we have to assume that in the solution of eq. (2.2.30) all the constant coefficients, 


analogous to the coefficient #,, in eq. (2.2.28), are identically equal to zero. An analo- 
gous procedure has to be applied to the neglected terms in eq. (2.2.17). 
Thus we have the result: 


Upstream [— co <a <0]: 


01 =0,+ Ey 0, U uj ar oe >) 2’ Es: 0, {exp [2 «—7 (1,4, y)]} sim (uez)- Ve; (2.2.31) 


me 


Downstream [0 <x <+ o]: 


012 = 02 +H, 0, U Tie 


yD Fin {exp [— 2 a+4 974 y]} sin (My 2)! uw ee), (2.2.32) 


u w 


2.3. Matching Conditions at the Disc 


For the remaining part of this analysis, we shall use the subscript 1 for x= —0, 
and the subscript 2 on x=-+0, these being the stations immediately upstream and 
downstream of the actuator disc. The choice of matching conditions follows those due 
to Yeh. The first matching condition at the dise is continuity. This is: 

(01+ @o1 + 1’) (U + toi + Uy’) = (02 + Goa + @2') (U + Ue t Ue’). (2.3.1) 


The second matching condition is that the spanwise velocity should undergo no change 
across the actuator disc, i. e., 


if 


Wy = War, (2.3.2) 
which implies that [see eqs. (2.2.6) and (2.2.7)]: 
Lewes Bons Uy’ ae Uy ; Oni == 05. (2.3.3) 


If the additional assumptions are made concerning the continuity of the velocity at 
the disc, i.e 
Ugy Uy =Upo + Us; (2.3.4) 
then 
Cr+ Cort C1 =O2+ Cont Ox» OF: Ogi t+ 01 =Oo2t Oe’, (2.3.5) 
if the assumption is made that ¢,=, which is justified by the fact that 0; (¢=1, 2) 


depend only upon the velocity U (very far upstream and far downstream, see eqs. 
(1.1.5a) and (1.1.5b). Let us assume another particular case, namely that m= jm = 0, 
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which leads to the result that w,’ and wy,’ are identically equal to zero and there are 
no conditions superimposed on the coefficients A,,, and B 


nm: 
wy’ =w,'=0. 2.3.6 
1 2 


This would imply that in this case the condition (2.3.3) is not valid. 


2.4. Cascade Characteristics 


Assume steady state conditions; a cascade is subjected to various relative inlet 
flow angles f,. One may consider the relative angle 6, as the independent variable, 
then on it there depend the relative exit flow angle £, and 
the loss of total pressure A p,. We shall consider at first 
the relative exit angle fy. as 

Consider only the kinematic velocity pattern; then 
following Yeh one can write: 


4 
U te 5 + vo ff ¢ 
tan (B+ A p)=—pr ee Cie LE 


But one may consider some additional elements like those Fig. 1. Velocity Pattern 
due to the downwash phenomena, etc. To allow for the 


variation of the ratio on the right hand side of eq. (2.4.1) one has to introduce a 
correction factor 6, i. e.: 


tan (B+ 6+A (6+ 6)) =[U tan (B+ 6)+0,+0'] [U+u+u'}. (a2) 
Apply a long division of polynominals to the right hand side of eq. (2.4.2): 
U tan (B+ 6)+u,+0': U+u,+w' =tan (6+ 6)+ (v)+’) U1 


—(uj+u') U-+tan (B+ 6)+ ..., (2.4.3) 
which may be written in the form: 
tan ((8 + 6)+ A (6+ 6)) = tan (6+ 6)+ A tan (6 + 6), (2.4.4) 
with: 
A tan (8 + 6)=(v)+’) U-!— (uy +u’) U- tan (6 + 6). (2.4.5) 
From: 
d/d (6 + 6) tan (B + 6) =sec? (6 + 0), (2.4.6) 
one gets: 
d (B + 6) =cos? (B + 6) d tan (f +6), (2.4.7) 
or ' 
A (8 + 6) = cos? (B + 6) A tan (8 + 0). (2.4.8) 
Using eq. (2.4.5) one has: 
A (B+ 6) = [eos? (8 + 6)] U- [(ep-+ 0") — (up +a’) tan (B+8)]. (2.4.9) 


Let us propose a relation between the exit flow angle, #2, and inlet flow angle, 6,, in 
the form: 
A (B.+ 6.)= WN’ A (8, + 9), (2.4.10) 


with 6, = 0. Using eq. (2.1.9) one obtains: 
(po + Vo’) — (og ty’) tan (By + d2)=N [(¥o1 + M1) — (Mor + M1) tan By], (2.4.11) 


with 


88 M. Z. v. Krzywoblocki: 


, N =N' cos? B,/cos? (B, + 9). (2.4.12) 
Using eq. (2.4.5) one can write: 

A tan (8, + 6,)=NA tan fj. (2.4.13) 

Using eq. (2.2.8) one can write eq. (2.4.11) in the form: 

Up tan (By+ 5g) —Vp2+(tan (By+ d2) + 4 Vem U = 

—N [uo, tan By — 91+ (tan By — 1% Ynm) Uy’). (2.4.14) 
This is the first fundamental characteristic relation of the cascade flow. As the next 
item let us consider the loss in the relative pressure across the cascade. The energy 


equation, originated from the first law of thermodynamics, has the following form 
(see eq. 2.1.12) at the point 1: 


4 (Oi + Vint Wh) +y (y—1)- pa Of —9 (81) —F (8) = 
=const. along a streamline, (2.4.15) 
with : 
On=0 uy tu: VaH=VitkhO 47421; Wu, (2.4.16) 
and similar expressions for U,2, etc. Applying the results of the shock theory to the 
present theory we assume that the constant in eq. (2.4.15) is the same on both sides 


of the actuator disc provided that one follows the same streamline. Apply the long 
division of polynomials to the expression: 


Pr1/ O11 = (Pr + Por t+ Pr')/(@1 + Cor + C1) = (P1 + Poi t Pr’) O1-* — 
(Piss Pora Pe) (Coa Oi OY ae (2.4.17) 


thus obtaining (after preserving only the first order terms in 1; 7’): 


Pirlo © (Pit Port P1 ) or" "Py (ax + Os) Oi: (2.4.18) 


Let us neglect in the expression (2.4.16) higher order terms like w,'?; wo, u,’; etc., thus 
obtaining from the first term in eq. (2.4.15): 


V242 0 (ayia) (Vath U2 Vy (gp 0) 22 BU oe). Cre 
Using the energy equation construct the difference: 


F1US 1 (VeteUyP I+ y y— 1) ple. $ (U2 Wy kU |) ee 
=U (ty + Uy) + (Vi +k U) (vp +04')— [U (Uo n+ Ua’) + (Voth U) (Vo2+e ) J+ 
+ ¥ (y— 1) [(o1 + Pr’) 01 — Pi (Cor + Or ) 01°] - 
—¥ (y—1)-* [(Poz+ Pa’) 02° — Do (Oo2+ 2’) 02°] — 
— [@ (81) +F (s1)] + [Q (82) +F (s9)]- (2.4.20) 


We assume that the left hand side of eq. (2.4.20) represents the required difference 
in the kinetic energy and the internal energy of the medium on both sides of the 
actuator disc. In the ideal conditions one suppose to get zero on the right hand side. 
In the real conditions one gets on this side an expression which represents the losses. 

Let us insert into eq. (2.1.5) the equality w’,,=v',,; into eq. (2.1.16) 0, =a 
and into eq. (2.1.7) w’,,=w',,; and w’,,=v’,,; let us multiply eq. (2.1.5) by da, 
(2.1.6) by dy, and (2.1.7) by dz, let us add these equations and integrate the result, 
thus obtaining the expression: 


Uy Voy Ey oe FH (2.4.21) 
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with a similar expression for the components with the subscript 2. Using eq. (2.4.21) 
in eq. (2.4.20) furnishes the following expression for the losses: : 


[U 1+ Vi %1+kU (Yori +21 )I—(U UpetVare+ hk U (Yo.+%5')] + 

+ y (y—1)™ [por 0: — Px (C01 + 01) Op et) pr pr — 

—¥ (y—1)™ [Poe 02 — Pe (Q02+ 02’) 02'] —(y— 1) pr’ Degas 

+ [Q> (8) —Q, (s)], (2.4.22) 


where the last term represents the total heat added to the medium during its passage 
through the actuator disc, the work of the external forces being excluded. The density 
magnitudes @, and g, must be known or calculated from the continuity equations using 
the velocity components U, V, and U, V,; the density magnitudes 0), and 095 and 
0, and e,’ were calculated in Section 2.2. Having these values, one can calculate 
approximately the values of the corresponding pressures, i. e.: 


Py! Poo = (01/@co)”s ~ Por! Poo = (001/@c0)*3 — (P1' Poo) = (011 O00)" (2.4.23) 


where p,, and o,, refer to any arbitrary standard conditions and the assumption is 
made that in the region 1 the flow isentropic. Similarly, in the region 2 the flow is 
assumed to be diabatic, i. e.: 


Pol P: = (@2/01)" exp [(S2—Sy) 67]; (Po2/Po1) = (@02/001)" * 
-exp [(S,—8,) 67]; pe’ py’ = (@2'/e1')” exp [(S2—S1) 7], (2.4.24) 
or 
P2|Poo= (02! Qc.)” exp [(S2—S..) 1]; Poa! Poo = (002/ @c0)” * 
“exp [(S,—S..) 7]; Pa" | Poo = (02'/Qc0)” exp [(S2—S..) *]. (2.4.24a) 
The entropy difference S,—S, refers to the heat addition (Q (s,)—@ (s,)) across the 


actuator disc. Following Yeh let us assume that one can express the losses across the 
actuator disc in the form: 


A & = losses/[4U? (1+ tan? B,)]=2 M (A B,), (2.4.25) 


where M is one-half of the slope of the @-function versus f, curve. Substituting 
eqs. (2.4.9), (2.4.22) and (2.2.8) into eq. (2.4.25) one obtains with (V,+k U)= 


= [tan (62+ 4)] 
| (t%o2+ + Voq tan (B2+ 02 ))— (491+ Vo1 tan By) —K 4 Y_m (Uy + Ug’) + 
+]UA=M [u,, tan B, —% 1+ (tan By —1 Yam) Uy’); (2.4.26) 
I=y (y—1)- [Boz 62* — Pe (Qo2 + @2') 027— Por Oi + Pi (C01 + @1') O17) + 
+ (y—1)** (p2' @2*— py’ @1") —(Qs (8) —@ (s)) - (2.4.27) 
Kags. (2.4.14) and (2.4.26) represent the cascade characteristics. Let us transform the 
expression J using eq. (2.3.5). Using the first eq. (2.4.24) with e,~ @,, we may expand 


the exponential function in a series and preserve only the first term, thus obtaining 
[with the use of eq. (2.3.5)]: 


Pr 01 — Pz Ox’ © Pr (01' — @2') © P1 (002— C01): 
Although p, +2, we can write in the first approximation: 


(D1 Co1— P2 Coz) + (Pi 01 — Po 02)» 
Dy (0o1— C02) + Pi (Co2— @01)* 9. (2.4.28) 
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Consider next the part of J equal to: 


(Po' @2'— Py’ Q1*) = (P2' — Pr’) Oi" (2.4.29) 
Using equations (2.4.23) and (2.4.24a) one gets: 
(Do! — Px’) = Poo Occ" {09" exp [(S2—S4) 1] — 01%}; (2.4.30) 


Por Dar Olo= 0%, EXP [(S2—S1) 654]; Por Pas ho = O13 Se—S.~8,—Si. (2.4.31) 
Following eq. (2.3.5) assume in the first approximation: , 
e2'% exp [(S,—S,) ¢1]— @1'7~ of, — of exp [(S2—S) 07]; (2.4.32) 


thus obtaining with the use of eqs. (2.4.31): 


(P2' — Py')* (Po1 — Poo). (2.4.33) 
This implies that with 0,~ 0.: 
Ix y (y—1)* (Po2— Por) @1' — (vy — 1)-* (Po2— Por) @1* — [Ge (8) —Q: (8) J = 
= 011 (Po2— Por) — [Qe (8) —Q (s)].- (2.4.34) 
One may attempt to justify this expression for 7. Assume in eq. (2.4.18): 
PrrlQa® Pili + Porleor + Pr |e’; (2.4.35) 


with the analogous expression for ,2/0,2. Let us use these expressions in the derivation 
of eq. (2.4.26) thus obtaining: 


Lm yy 1) (Pog Ore Pa 03 Pol Oi Ps a 
— px’ 021 — [Qe (8) —Q, (s)], (2.4.36) 


or with @2% 023 01% 015 02* 01: 


Ixy (y—1)* (Po2— rede = (Dal = O19) C1 
=o(y a Pit) — [2 (s) —Q, (s)]. (2.4.37) 
Use of the result (2.4.33) in eq. hohe) furnishes exactly eq. (2.4.34). 


One may provide another proof of the validity of eq. (2.4.33). Namely, let us 
begin with the momentum theorem in the horizontal direction when no forces are acting 
in that direction: 


Pat ou Ui = pz + Q2U io, (2.4.38) 
or 


(Pi+ Por + Pi) + (01+ O01 + O1') (U + U9, +4’? = (0, + Poot Pe)+... (2.4.39) 


Neglecting second order terms like u61, u,', Uo, U1’, 01’ Uo, ONE gets: 01! Uy’, O91 Uor Cor 
U,', etc., 
Pit Q, [U? +2 U (ug, + %y')]+ (001 + 01’) 02 © 
~ Pr2+ Oo [U2+ 2 U (Uye-+ Ue’)]+ (092+ 09’) U?. (2.4.40) 
Taking into account oS (2.3.4), (2.3.5), the condition 9,~* 02, and the assumption 


that in general, p);, p;’ (t=1, 2) are small in comparison with p;(t=1, 2) which 
implies that p,~ p., or 


Port Pr © Poot Pe’, (2.4.41) 
which is equivalent to eq. (2.4.33). 
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3. Outline of the Solution of the Problem 
1. General Philosophy 


Assume that we know all the functions in the region 1 (uw 1, %1, %, By, ete.) except 
U;’, V;' and w,'; in the region 2 we know a», By. Using eqs. (2.4.23) and (2.4.24) we can 
always express the density @ in terms of the pressure p and vice versa. We shall also 
use conservation of energy across the actuator disc with the assumption from the 
shock-theory that the constant in eq. (2.4.15) is the same on both sides of the actuator 
disc on one and the same streamline. Assume that the following values are given 
on the side 1: 


Mh Opa (24223) 55 teas 057 (22:23)3. Po (24.23)% 
(2.2.15 


a1 ( 
tatigee Gy, V7 UG. 04 (2.2. 15):, tan 6, 0, ey Ba 


on the side 2: 
Us Q2=013 Q (52) —@ (51) =Qs (8) —Q (8); Pe (2.4.24); 
F (s,)—F (s,)=F,—f,; M;-N; tan (6:4 6,)=0,, (3.1.2) 
and the following 8 unknown values are sought: 
On the side 1: 
: Mee, Og Wg (3.1.3) 
on the side 2: 


oo (and 095 (2.2.24) and py. (2.4.24)); 
Vas, (consequently, talire,—=d,);. Us U5 se ie (3.1.4) 


We do not mention the quantities 0,’ and 0,’ [and consequently p,'; p.'; (2.4.24) 
and (2.4.24a)| since using eq. (2.3.5) one may attempt to represent the quantities 
0;'; Qo’; in terms of 09, and 9y,. To our disposal in the case of m +0, umn +0, we have 
a system of 7 equations: 


(2.2.8) — 2 eqs.; (2.3.2); (2.3.3); (2.3.4); (2.4.14); (2.4.26). 
In the case m =u», =90 (which implies that w,’=w,'’ =0), we have the following 
system of 5 equations: 


(2.2.8) — 2 eqs.; (2.3.4); (2.4.14); (2.4.26), m,6 unknowns 
Lies Se MDa Gs (OMA Ogg BIG SRN Dios, ha SSA 


3.2. Particular Case of Spanwise Distortion 
This case is characterized by: 
Mee 0 eee eO. (3.2.1) 
Eqs. (2.2.8), (2.3.2), (2.3.3) and (2.3.4) give, respectively: 
Chea 056 Os 5 hy hg 5 tg ps — 2 ty S20 (3.2.2) 
Eq. (2.4.14) gives: 
(6, —N 6,) Ug — 2 Ugg =(N U1 — be) Up — 2 NU 1. (3.2.3) 
With the use of equations (3.2.3) and: 
(Poal'P2) = (00 2/02)” xP [(So2—S2) 6511; Soe © Se: (3.2.4) 
(Po1/P1) = (@01/01)"; CtC.; or Uos= Le Oy U; ete., (3.2.5) 
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eq. (2.4.26) gives after all the manipulations are performed: 
[(bg—N by) bg— Mb, +2] wh, wis? + 
A [— Bath — 2 ad! tog Op 
— 2 U-1 (Q, (s) —Q, (s)) u%, —(N Uo — 02) Utt* be + 
+2N uk, 091 b,.—M utf* 6, +2 M uf, voi] ufo + 
+ 2 HY UY-1 {exp [(So2—S2) al} Po Qi Uh, =9; 
Soo Se; Soe—So~ 09. (3.2.6) 


To solve eq. (3.2.6) one may propose an iteration process. With y=1, one obtains 
‘a second degree algebraic equation which readily can be solved, giving result to be 
denoted by the symbol woo (9). Next, one can apply an iteration process of the kind: 


[...] wei) =[---] Wort) + free term, n=1, 2,3... (3.2.7) 


The iteration process may be carried on up to the required degree of approximation. 
Eas. (3.2.3) and (3.2.6) determine two unknowns wy». and vy5. The magnitude p, may be 
substituted by p, using the equation: p,/p,= (0s/01)” exp [(S2— 4) ¢,1]- 


3.3. Circumferential Distortion 


In this case one has: m=pm=0; %am=1. This is the second case discussed in 
Section 3.1. Eqs. (2.2.8), (2.2.15a) and (2.3.4) give: 
0, 20 yp Ue = SS ee es (3.3:1) 
Vp o = Ups tan a. (3.3.2) 
iq. (2.4.14) gives: 
uy’ = [(N k—b,—14) Uy + (A. +14) Upe] * 
(6, — Nb.) 0 (3.3.3) 
and eq. (2.4.26) furnishes: 
US, (Cott Do) user? —{usit? (Ay +4 Bo) + By’ U™ pyloy + 
+ U-* ub, [Qe (8) —Q: (8) ]}uss 
+ HY UY ut, {exp [(So2—S2) ¢1]} Ps 11 = 0. (3.3.4) 
Excluding the imaginary terms in eq. (3.3.4), one obtains an equation for uy,. The 
magnitude wu,’ may then be found from eq. (3.3.3). The coefficients are: 
Ag= (1+ by + bk) (b,—N b,) + M b, (N k—b,)+ (kh —k N—M); 
By=(1+a,6,+Mk)1+N)+kN (2k—b,—M)—kb, + M (b,—6,); 
Co= (1+ a2 6.) (6, —N b,)-+k (1— N)— M (1+, b,); ; 
Dy=1+a2+kh?+ M (a,—b,)+ N (l+a,b,—k),). 
(Received September 29. 1959) 
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(Staatliche Forschungsanstalt fir Warmetechnik, Praha) 


Kinflu8 der elastischen Fundamentlagerung auf die Rotorstabilitat 


bei Beriicksichtigung der inneren und 4uferen Dampfung 
Von 
A. Tondl, Praha 


Mit 4 Textabbildungen 


Zusammenfassung. In der Arbeit wird der Einflu8 des gleichzeitigen Einwirkens innerer Werk- 
stoff- und 4uBerer Dampfung eines Turbinenliufers wie auch der Dampfung im elastisch gelagerten 
Fundament an der Grenze der stabilen Bewegung dieses Systems analysiert. Weiter wird der Grenz- 
wert der Rotordrehzahl bestimmt, oberhalb dessen die von Unwuchten erregten Rotorschwin- 
gungen unstabil werden. 


I. Einleitung 


Der KinfluB der inneren Dimpfung (Materialdimpfung und Dampfung, verursacht 
durch die Reibung miteinander verbundener Teile bei der Deformation eines Rotors, 
z. B. Reibung zwischen Welle und Nabe der Scheibe) wurde bereits von einer Reihe 
von Autoren!-* analysiert. Fast alle Arbeiten beschaftigen sich mit dem Falle, 
in welchem der Rotor auf ideal starren Stiitzen gelagert ist. Eine Ausnahme machen 
zwei Losungen, welche in den Arbeiten 7 und 8 enthalten sind. In der ersten Arbeit 
wird der Kinflu®B masseloser elastischer Stiitzen und in der zweiten der Einflu8 der 
elastischen Lagerung des Fundamentes beriicksichtigt. In dieser zweiten Arbeit® 
nimmt der Verfasser F. M. Dimentberg die innere Dampfung als lineare Funktion 
der relativen Geschwindigkeit an. Die Losung fiihrt mit Riicksicht auf die Linearitat 
der Differentialgleichungen zu den Hurwitzschen Stabilitatsbedingungen. Da die 
praktischen Schwierigkeiten bei der Berechnung der Hurwitzschen Determinanten fiir 
eine allgemeine Losung der Aufgabe zu groB sind, beschrankte sich Dimentberg blo8 
auf einen speziellen Fall, wenn die Massen des Fundamentes und der Rotorscheibe, die 
Koeffizienten der inneren Dampfung und der Fundamentdampfung und die Eigen- 
frequenz des Rotors und des Fundamentes iibereinstimmen. Auf diese Weise konnte er 
jedoch nicht den EinfluB der Abstimmung des Fundamentes, den Einflu8 des Verhalt- 
nisses der Rotoren- und Fundamentmassen usw., analysieren. 

In der vorliegenden Arbeit werden wir die Analyse eines ahnlichen einfachen 
Systems durchfiihren, wobei wir uns jedoch nicht auf die Ubereinstimmung der Para- 
meter des Systems beschranken. 


1A. L. Kimball: Internal Friction Theory of Shaft Whirling. General Electric Review, 
April, 1924. 

2D. Robertson: Hysteretic Influence of the Whirling of Rotors. Proc. of the Inst. of Mech. 
Eng., vol. 131, 1935. 

3 F.M. Dimentberg: Ob ustojéivosti gibkogo vala s neuravnoveSenym diskom pri dejstvii 
ynutrennego i vnesnego trenij (Uber die Stabilitit einer elastischen Welle mit nicht ausgewuchteter 
Scheibe bei Einwirken innerer und iuRerer Dampfung). Izvestija A N SSSR, ot. tech.n., No. 10, 1954. 

4M. I. Cajevskij: Issledovanije ustojéivosti vraSéenija gibkich valov za kritiéeskoj skorostju 
(Forschung auf dem Gebiet der Stabilitiat elastischer Wellen bei tiberkritischen Drehzahlen). 
Izvestija A N SSSR, ot. tech. n., No. 9, 1955. 

5 E. Downham: Theory of Shatt Whirling. The Engineer, October 11, 1957. 

6 A. Tondl: Der Einflu8 der Werkstoffdimpfung auf die Stabilitat der Rotoren. Acta tech- 
nica, No. 3, 1959. 

7A. Tondl: Nékteré otazky kmitani a stability odpruzenych rotort (Einige Fragen der Schwin- 
gungen und Stabilitit elastisch gelagerter Rotoren). Rozpravy CSAV, tada TV, Ne. 5, 1955. 

8. M. Dimentberg: Popereényje kolebanija vraséajuséegosja vala pri nalicii soprotivlenija 
trenija (Biegeschwingungen einer rotierenden Welle bei Einwirkung der Dampfung). Popereényje 
kolebanija i kritiéeskije skorosti, sbornik pjervyj, Izdat. A N SSSR, 1951. 
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Ausfiihrlichere Untersuchungen iiber den Hinflu8 der inneren Dampfung zeigten’, 
daB8 — sofern nur die Stabilitat der Bewegung untersucht wird und ein even- 
tuelles Aufkommen von selbsterregten Schwingungen nicht von Interesse ist — die 
Resultate der Losung mit nichtlinearen Funktionen der inneren Dampfung qualitativ 
mit der linearen Approximation iibereinstimmen. Fiir unsere qualitative Analyse wer- 
den wir daher die innere Dampfung als lineare Funktion ansetzen. Auch die aubere 
Dampfung werden wir mittels linearer Funktionen approximieren. 

Wie bereits erwahnt wurde, ist die allgemeine exakte Bestimmung der Stabilitat 
linearer Systeme mit mehreren Freiheitsgraden, welche zu den Hurwitzschen Bedin- 
gungen fiihrt, praktisch undurchfiihrbar. Daher werden wir unsere Aufgabe mittels 
einer Naherungsmethode lésen. Diese Methode erlaubt uns, auf bequemere Weise die 
Stabilitétsbedingungen mit beliebiger Genauigkeit festzustellen und ermdglicht, die 
Analyse aller oben angefiihrten Einfliisse durchzufiihren. 


II. Die Methode 


Es sei gegeben ein System linearer homogener Differentialgleichungen mit kon- 
stanten Koeffizienten 


an +X st weed (Pas Lk + Csk Xx) (s= 1, 2, 3, enete:s n), (1) 


wo ¢ ein kleiner Parameter und a,x, b;4, ¢,, Konstanten sind. Das im weiteren beschrie- 
bene Verfahren ist vor allem dann giinstig, wenn das gektirzte System (1) 


n 


° Ls Se Ask LE = 0 (2) 
k=1 


in mehrere unabhangige Systeme zerfallt, so da die Eigenfrequenzen des gekiirzten 
Systems Q, (k= 1, 2,..., 2) leicht zu bestimmen sind. Man nehme an, die Wurzeln der 
charakteristischen Gleichung 


det | Ask + Osn A | a) (6s.= 0, S = Taye One a5 1) ‘ (3) 


waren alle rein imagindr und voneinander verschieden: &% = + iQ; i= \—1; 
Aree ApeihS Hylt) 5 1g" =O. 


Wir transformieren das System (1) derart, daB das gekiirzte System (2) in normale 
Koordinaten iibergeht und erhalten so ein System 


Ys— AP Ys = 5 (Dek Ye Qe Yr) (SAU, 2,3, 4 0). (4) 
Mit Hilfe der Substitution . 
Ys As Ys = Us, Ys — As Ys = Us 4 &(D) 
fiihren wir das System (4) in die Form 
ti = Ate +5 eP, 


‘ 1 
Un. — ApOe- Ss 6s _ (Saas 2ieden- ce, 


iiber, wobei 
n 


‘a 1 
dhe “ [Poe (us + vy) + 4 Usk (Ur — Un)|. 
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Da das System (6) linear ist, konnen wir eine Lésung in der Form 
Mas Lie pss Vic ett (7) 


suchen, wo U,, V, von der Zeit unabhangige Konstanten sind, die jedoch den Parameter 
e enthalten. Der charakteristische Exponent » ist gleichfalls eine Funktion des Para- 
meters ¢«. Da jedoch » = 4 gilt fiir ¢ = 0, kann y in der Form 


y=A+t €a(e) (8) 
geschrieben werden, wo 2 die Wurzel der charakteristischen Gleichung (3), z. B. 4; (j= 1, 
2,...,%), und « eine gewisse Funktion von « sind, wobei bei « vorausgesetzt werden 


kann, daf§ eine Entwicklung in den Potenzen von ¢ existiert. Man suche daher die 
Lésung in der Form 


u;—=AUg +e Us +e Ung baw. Jem beta t that «dt, 
, 9 
Vs = (Vso te Vireo? Vaat...) Ut eem tent .I (9) 


Durch die Substitution von (9) in (6) und durch den Vergleich der Koeffizienten glei- 
cher Potenzen von « erhalten wir fiir die erste Annaherung 


U0 (A, — 4;) = 0, 
V;,0 (As +4) =0 UR econ Ee a A 


eeeeiee eae Wan Op, = 0) (haa 12,59, 4/03 ke), Vio Ors = 1,2, 330 5 st). 
Der zweiten Approximation entnehmen wir 


(10) 


1 1 
%,1 = g (Dis FG, Ui 


Uri=35 Ga te A = VAD oe on 105 ela (11) 


} Se 24 Aj + As 
Aus der dritten Approximation kann bestimmt werden 


1 ij 1 1 Sapte 
OG as [pie Re =a (a ny 4:) b Uk Fa ye (a 5 4) (12) 
cit 


Falls wir uns mit der Annaherung ‘ 
p= Ay + & 1 + 8 o,2 


begniigen, ist es nicht notwendig, weitere Werte U,.,,, V.,, und %,, zu berechnen. Die 
Stabilitat der Lésung ist durch die Bedingungen gegeben, da} die realen Teile aller 
Cag l, 2. en) negativ sind : 


Rete) = UL tes, 2.0): (13) 


wobei «; = a1 + € 4,2 +... mit beliebiger Genauigkeit bestimmt werden kann. 


III. Die eigentliche Lésung der Aufgabe 


Es sei ein einfacher Rotor mit einer in der Mitte der Welle aufgesetzten Scheibe 
gegeben, dessen Masse gegeniiber der Masse der Scheibe m vernachlassigt werden kann. 
Das Fundament mit einer Masse mp sei zur Ebene der Scheibe vollkommen symmetrisch 
aufgebaut, und zwar sowohl in bezug auf die Verteilung der Masse als auch in bezug auf 


Wf 
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die Elastizitat der Fundamentlagerung. Falls angenommen wird, daB die Scheibe auf 
der Welle nicht schrag sitzt, wird die Bewegung der Scheibe und der Fundamentmasse 
in einer Ebene verlaufen. Die Steifigkeit der Fundamentlagerung sei sowohl in verti- 
kaler wie auch in waagerechter Richtung ¢,. Die Steifigkeit der Welle sei c. Weiter sei 
angenommen, der Rotor drehe sich mit konstanter Winkelgeschwindigkeit o. 

Wir bezeichnen nun die Schwerpunktkoordinaten der Scheibe mit x, y und die des 
Fundamentes mit 2, y). Es sei zweierlei 4uBere Dampfung betrachtet: die Rotor- 
dampfung, proportional den Geschwindigkeiten « — 2), y — y, mit einem Koeffizienten &, 
und die Fundamentdimpfung, proportional den Geschwindigkeiten x), y, mit einem 
Koeffizienten k,. Die innere Dampfung, wie bereits gesagt, werden wir mit Hilfe einer 
linearen Funktion der relativen Geschwindigkeiten im rotierenden Koordinatensystem 
approximieren (das Koordinatensystem rotiert mit dem Rotor). Im festen Koordinaten- 
system ist daher die Dampfung den Geschwindigkeiten 7+ oy und y—wx mit dem 
Koeffizienten h proportional. 

Auf Grund dieser Erwaigungen kénnen die Bewegungsgleichungen leicht aufge- 
stellt werden. Da die Gleichungen linear sind, kann fiir die Berechnung der Stabilitat 
das homogene System betrachtet werden: 


m £4 ¢(%— Xo) + (k-+h) (%— xo) + wh (y— Yo) = 9, 
my +e (Y¥—Yo) + (k+h) (y¥— Yo) — wh (x — ao) = 0, 


: a ee (14) 
Mo Xo + C (Ly— 2X) + Co Ly + ky Z + (K+h) (%o— Xo) + Wh (Yo—y) = 9, 
Mo Yo + € (Yo— Y) +60 Yo ko Yo (+4) (Yo— y) — wh (Xp — 4) = 0. 
Falls wir die Bezeichnung «+7 y = z einfiihren, erhalten wir nach Umformung 
2+ — (z—2z) = — e[(x + 0) (2%) iw 6 (229) ], aes 
15 
Zo+ pe (20-2) + to = — 6 [mo Zo + (4 + 44) (2 — 2) — 1 Oy (%—2)], 
wo 
Oo yk De sae _>~h ees ee 
Oe 3 A haga E ae EO, mo” xo ™mo 
Mit Hilfe der Substitutionen 
Z—= Uy + Ug, % = AU, + 6 Up, (16) 
wo 
c ¢ 
mo mo 
a= ; = 
e + Co O,2 e+ 0 = Ore : 
mo Mo a 
Fee lh G c+ co / ¢ \2 c+ ¢0\2 2¢(e€— eo) 
21,2 4 = Mo + \/(5) | ( mo ) mmo 
ist, iberfiihren wir das System (15) auf die Form 
uy — Ay? uy = € (Pay ty + Pye Ue + Qi Uy + Qo Ue), 
Uy — Ag? Uy = € (Poy Uy + Poe Ue + Qo Uy + Qos Ue), (17) 


wo 
Ax,2 == 0) Qy,2, 


Pyy= pz {@—1) [0 (8-+%) + dy tm] + am, 
Prepay (01) 1B +0) + 6, bond + b ml, 


Pu= pz {(1—a) [a (6-+%) + 8, 4%] —a xg}, 
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1 
Pro= p—G {(L—}) [a (6 +x) + 61+] — bx}, 
Qu=p—gio(db+6,)(1—a), Qe= 
Qu = p—qio(da+ d,)(a—), Caer 1 @ (da+ 64) (6—1). 


Durch Anwendung der oben angefiihrten Methode erhalten wir 


f 


para COB 401 3) 


Re (a,1) = 5g \(@—1) [6 (6+%)+ d+%]+4 not G- ~ (1 —a) (6 6+ s)}. 

Re (1,2) = S (18) 
] { 7) 

Re (2,1) = p—q (1-9) [a (6 +2) + 8, +4] — b m+ 2 (6-1) (a d+ 6,)}, 

Re (a 2) = 


Falls wir uns mit der Annaéherung 


y=A+ea,e? a, 
begniigen, erhalten wir folgende Bedingungen fiir die Bewegungsstabilitit: 


o<O(1+F54+2K i) 


(19) 
o<2,(14+54+3 K,), 
wo 
: a 1 ; b 1 
Ky =p rap me K,=574 Ta tie ist. 
m 

Ks kann nachgewiesen werden, dai folgende Gleichung gilt 

aps ay (20) 


Daher wird 
K,= (25). Ke=(55)- 


Die Stabilitatsgrenze ist gegeniiber den Kigenfrequenzen 2, und 2, erhéht worden, 
einesteils durch das Einwirken der auBeren Dampfung des eigentlichen Rotors und 
anderseits durch die Dampfung der Fundamentmasse. Der Anteil der auBeren Rotor- 
dampfung an der StabilitatsgrenzenerhOhung hangt direkt vom Verhaltnis (k/h) der 
Koeffizienten der au8eren und inneren Dampfung ab. Der Anteil der Dampfung des 
Fundamentes hingt einesteils vom Verhaltnis (ko/h) der 4uBeren und inneren Dampfung 
des Fundamentes, anderenteils von den Ausdriicken AK, und K, ab, die ihrerseits von der 


Abstimmung des Fundamentes und vom Verhaltnis der Massen aL abhangig sind. 


Da auch 2, und 2, von der Abstimmung des Fundamentes abhangig sind, wird es 
fiir die Schatzung der Vorteilhaftigkeit der Abstimmung des Fundamentes giinstiger 
sein, die Grenzwinkelgeschwindigkeit des Rotors als dimensionslose Grofe, namlich 
als Quotient der Winkelgeschwindigkeit des Rotors und der Eigenfrequenz des Rotors 
auf absolut starren Stiitzen auszudriicken: 


(0) Qy k ko a 2 Qe k ko b 2 
> PS ears f eae (5)'|, ae erg eg) |} (21) 


wobei @, = ee die Kigenfrequenz des Rotors auf absolut starren Stiitzen ist. In Abb. 1 


ist das Grenzverhaltnis 7 in Abhangigkeit vom Abstimmungsverhaltnis __, dargestellt, 


wobei wy = en die Eigenfrequenz des Fundamentes ist, und zwar fiir = =z Ly 
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| ieremrt bP =aee — 1, Wenn wir das Resultat mit dem Fall der absolut starren Lagerung 


(7 = 1) vergleichen, dann sehen wir, da® durch das Einwirken der elastischen Lagerung 
und der Fundamentdimpfung das Grenzverhaltnis stets groBer als 1 ist. In Abb. 2 


3 
7 
2 
4 
(6) 74 
a 1 Sick Be 
@, ay 
Abb. 1 Abb. 2 
of See ko/mo = 
ist itr Lk OL dargestellt. Aus dem Verlauf ersehen wir, dai das 


kleinere Verhaltnis ~ sich giinstiger auswirkt. 


«4 = pe, Mo k > ‘ko/mo : awe 
In Abb. 3 ist 7 fir —- = 10, 5 = 2, eat und in Abb. 4 fiir —) 


no 


k 
Liv Saeae 
==" igcat= 


Abb. 3 Abb. 4 


~~ = 5 dargesteilt. Aus dem Vergleich mit dem Fall eines auf absolut starren Stiitzen 
gelagerten Rotors, wobei 7 = 3 ist, ersieht man den giinstigen EinfluB elastisch gelager- 
ter Fundamente. 

In gr6Btem Ausmafe wirkt sich die Fundamentdimpfung dann aus, wenn die Kigen- 
frequenz des Fundamentes mit der Eigenfrequenz des Rotors iibereinstimmt. 
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Falls jedoch diese Eigenfrequenzen nicht iibereinstimmen, dann ist es giinstiger, 
wenn die Higenfrequenz des Fundamentes niedriger ist als die Eigenfrequenz des Rotors. 

Die Ergebnisse kénnen in folgende Punkte zusammengefaft werden: 

1. Die Grenze der Rotorstabilitat ist bei elastischer Lagerung des Fundamentes 
immer hoher als bei starrer Lagerung (bei Beibehaltung des Rotordimpfungskoeffi- 
zienten). 

2. Die Hebung der Grenze ist desto starker, je gréBer die Dampfung des Fundamentes 
ist. Die Kinwirkung der Fundamentdampfung macht sich bei niedriger Abstimmung 
des Fundamentes (Kigenfrequenz des Fundamentes niedriger als Eigenfrequenz des 
Rotors) mehr bemerkbar als bei hoher Abstimmung oder bei kleinerem Verhiltnis der 
Fundament- und Rotormassen (jedoch bei der Annahme gleichbleibender Fundament- 
dampfung). Bei der Konstruktion eines leichten und elastischen Fundamentes ist daher 
auf geniigende Dampfung zu achten, insbesondere bei Stahlfundamenten. Geniigende 
Dampfung kann durch Ausniitzung der Konstruktionsdampfung in den Verbindungs- 
stellen einzelner Teile erreicht werden. 

Diese Dampfung wirkt sich auch giinstig auf die Herabsetzung der Schwingungs- 
amplituden der Maschine aus. 

(Hingegangen am 4. Februar 1960) 


Der Einflu8 des Kriechens auf die Tragfihigkeit 
von Stahlbetonsdulen* 
Von 
Natalija Naerlovié-Veljkovité, Beograd 
Mit 26 Textabbildungen 


Die Frage des Einflusses des Kriechens auf die Tragfaihigkeit von gedriickten Staben 
wurde in den letzten Jahren in der technischen Literatur vielfach behandelt, und zwar 
im Hinblick auf Flugzeugkonstruktionen unter hoheren Temperaturen. Dieselbe Frage 
tritt aber auch bei Stahlbetonsaiulen in der Baupraxis auf. In beiden Fallen handelt 
es sich um die gleiche Eigenschaft des Materials, dal} namlich dessen Verformungen 
unter konstanter Last mit der Zeit zunehmen und das Spannungs-Dehnungs-Gesetz 
nicht mehr linear ist. Zur Beschreibung des Verformungsvorganges eines solchen 
Materials wird die Zeit als neue Variable in das Stabilitatsproblem von auf Druck 
beanspruchten Konstruktionsteilen eingefiihrt'. Damit erweist es sich als notwendig, 
den Begriff der ,,Stabilitat‘‘ zu erweitern. Unter dem in der amerikanischen Literatur 
iiblichen Terminus ,,Creep Buckling“ versteht man nicht mehr eine Eigenwertaufgabe 
im klassischen Sinne, sondern die Frage nach der ,,Lebensdauer‘‘ der Konstruktionen: 
wie lange wird es dauern, bis die Querverschiebungen eines mit baupraktisch unver- 
meidlicher Anfangskriimmung ausgefiihrten, axial belasteten Stabes so groB geworden 
sind, daB die Tragfaihigkeit des Materials erschopft ist ? 


_ * Die Anregung zu dieser Arbeit stammt von Herrn Professor H. Parkus. Fir die vielen wert- 
vollen Hinweise und Diskussionen wiihrend der Ausfiihrung dieser Arbeit an seinem Institut an 
der Technischen Hochschule Wien méchte ich Herrn Professor Parkus meinen herzlichsten Dank 
zum Ausdruck bringen. 

Die numerischen Rechnungen wurden auf einer elektronischen Rechenmaschine, Type IBM 650, 
im Mathematischen Labor des III. Institutes fiir Mathematik der Technischen Hochschule Wien 
durchgefiihrt. Dem Vorstand des Institutes, Herrn Professor Dr. R. Inzinger, mochte ich auch 
hier meinen besten Dank aussprechen. 

1B. Fraeijs de Veubeke: Creep buckling. In: High Temperature effects in aircraft struc- 
tures, ed.: N. J. Hoff, Pergamon Press, 1958. 
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Die Losung dieser Aufgabe erfordert Untersuchungen, die in zwei Richtungen 
zu fithren sind. Erstens ist fiir ein gewisses reelles, experimentell untersuchtes Material 
das zugeordnete ideelle, durch ein entsprechendes Spannungs-Dehnungs-Gesetz 
charakterisierte Material zu finden, das sich dem reellen Material méglichst gut anpaBt. 
Zweitens ist bei gegebenem Spannungs-Dehnungs-Gesetz, geometrischen Daten und 
Belastung des Stabes das Problem analytisch zu lésen. 

Die Beziehung zwischen Spannungen und Verformungen soll mit experimentellen 
Ergebnissen iibereinstimmen und soll den Bedingungen der mathematischen Theorie 
des isotropen unelastischen Kontinuums geniigen. In dieser Hinsicht unterliegt die 
Freiheit der analytischen Formulierung von experimentellen Ergebnissen mancher 
Beschrankung?. 

Wegen des ziemlich umfangreichen mathematischen Apparates, der schon im Falle 
des beidseitig gelenkig gelagerten Stabes zur Behandlung von Aufgaben dieser Art 
notig ist, blieb die Anzahl der bis jetzt gelésten Aufgaben klein, obwohl in neuester 
Zeit auch Lésungen auf dem Gebiete der Flichentragwerke erschienen sind**». 
An Lésungen von ,,Creep-buckling“‘-Problemen gelenkig gelagerter Stabe ist eine Reihe 
von Arbeiten bekannt, welche sich auf den idealisierten H-Querschnitt beziehen*. Uber 
den Rechteckquerschnitt scheint es nur vier Arbeiten zu geben. 

Charles Libove® diirfte der erste gewesen sein, der sich mit einem Stab mit 
rechteckigem Querschnitt befaBt hat. Die Form des gewahlten Spannungs-Dehnungs- 
Gesetzes ist in Ubereinstimmung mit dem Verhalten von Aluminiumpriifstaben 
gewahlt und entspricht einem visko-elastischen Material. Fiir eine approximative 
Lésung der Aufgabe beniitzt Libove Interpolation (collocation-method) zur Elimi- 
nierung der Veranderlichen in Richtung der Stabachse und verwendet weiter numerische 
Verfahren. Eine Kritik seiner Arbeit findet sich in dem Buch von N. J. Hoff**. 

Den Fall des Stabes mit rechteckigem Querschnitt behandelt weiters Sharad 
A. Patel in zwei Arbeiten. In der spateren’ wird die Lebensdauer fiir visko-elastisches 
Material mit mehreren Varianten des Spannungs-Dehnungs-Gesetzes berechnet. Die 
erhaltenen Ergebnisse fiir den idealisierten H-Querschnitt werden mit den Ergeb- 
nissen fiir den Rechteckquerschnitt verglichen. Wie bei Libove, wird auch hier 
Interpolation in Richtung der Stabachse und numerische Integration in der Quer- 
richtung und nach der Zeit verwendet. 

Die Arbeit von Bleich und Dillon® bezieht sich wieder auf ein visko-elastisches 
Material mit dem Ziel der approximativen Bestimmung der oberen und unteren Grenze 
der Lebensdauer. 

Die Anfangskriimmung ist in allen erwihnten Arbeiten in Form einer Sinushalbwelle 
angenommen. Kine Diskussion des Einflusses der hdheren Harmonischen fiir den Fall 


2 A. M. Freudenthal and H. Geiringer: The mathematical theory of the inelastic con- 
tinuum. Handbuch der Physik, (S. Fliigge), Band VI, Elastizitat und Plastizitat (1958), S. 257. 

3 F. W. French and A. Sh. Patel: Creep buckling of cylindrical ae subjected to uniform 
axial compression. PIBAL Report No. 489, May 1959. 

4 A. Sh. Patel, K. A. V. Pandalai and B. Venkatraman: a sunas analysis of thin- 
walled structures. PIBAL Report No. 497, July 1959. 

° L. Sanders and Harvey Jr.: A variational theorem for creep with application to plates 
and columns. National Advisory Committee for Aeronautics, Report 1342 (1958). 

* s. FuBnote 1. 

®° Ch. Libove: Creep-buckling analysis of rectangular section columns. NACA, June 1953. 
Technical Note 2956. 

** g, FuBnote }. 

7 A. Sh. Patel and B. Venkatraman: Creep behavior of columns. PIBAL Report No. 422, 
May 1959. 

* H. H. Bleich and O. W. Dillon: Nonlinear creep deformations of columns of rectangular 
cross-section. Journal of Appl. Mech., December 1959, S. 517. 
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des Stabes mit idealisiertem H-Querschnitt findet sich in der erwahnten Arbeit von 
Sharad A, Patel*. 

In der vorliegenden Arbeit wird versucht, die Lebensdauer eines gedriickten Stabes 
aus Stahlbeton mit rechteckigem Querschnitt, gelenkig gelagerten Enden (Abb. 1) 
und mit der Anfangsdurchbiegung 


Wy =%, sina é (a) 


zu bestimmen. Das ideelle Material ist ein quasi-lineares elasto-viskoses Material mit 
dem Spannungs-Dehnungs-Gesetz 


é=H(I—xo) +A WZ, (b) 


wobei der veranderliche Koeffizient eine lineare Funktion der Spannungsinvariante ist. 
Bei der Formulierung des Gesetzes wurde von einer Spannungs-Dehnungs-Linie 
ausgegangen, wie sie Ergebnissen von Versuchen an Stahlbeton entspricht, die an der 
Technischen Hochschule Hannover ausgefiihrt wurden’. Das dort gegebene Gesetz 
wird folgenderweise geschrieben 
oO x% 
é= EB (1 aa O O ), 


wo x ein Festwert ist. Differentiation nach der Zeit liefert den ersten Teil des Aus- 
drucks (b). Der viskose Teil des Gesetzes (b) soll den Kriechvorgang des Betons be- 
schreiben. Die Form der Kriechfunktion 4 (¢) wird nach Pucher" gemaB experimentellen 
Daten angenommen. Sie ist eine Exponentialfunktion der Zeit von der Form 
AG) =A.eF* 

wo A und B Materialkonstanten sind. 

Das Gesetz (b) ist nur fiir Betondruckspannungen und nur im Belastungsbereich, 
o > 0, giiltig. Es wird aber zunachst auch dem Entlastungsfall, ¢ < 0, zugrunde gelegt 
und auch im Zugspannungsbereich verwendet. 

Die Saule wird langsbewehrt und vorgespannt angenommen. Es wird vorausgesetzt, 
daB fiir den Bewehrungsstahl das Hookesche Gesetz gilt. 

Der Kiirze halber soll im folgenden statt der Bezeichnung ,,Ideelles Material“ stets 
der Ausdruck ,,Beton‘‘ verwendet werden. 


Bezeichnungen 
x — Koordinate in Richtung der Stabachse 
z — Koordinate quer zur Stabachse 
1 — Lange des Stabes 
2b — Hohe des Stabquerschnitts in z-Richtung 
' ¢ — Breite des ,Stabquerschnittes 
F — Querschnittsflache 
F, — gesamte Fliche der Langsbewehrung 
hm = — Bewehrungsanteil 
e — Abstand des Bewehrungsschwerpunktes von der Querschnittsachse 
f= + — dimensionslose Koordinate in «#-Richtung 
ee = — dimensionslose Koordinate in z-Richtung 


* gs. FuBnote’. 

°K. Gaede: Knicken von Stahlbetonsiulen unter Kurz- und Langzeitbelastung. Deutscher 
Ausschu8 fiir Stahlbeton, Heft 129, S. 54. 

10 A, Pucher: Lehrbuch des Stahlbetonbaues, 2. Aufl., Wien: Springer- Verlag, 1953, S. 20—21. 
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wy — Durchbiegung der Stabachse in z-Richtung fiir den unbelasteten Stab 


w — gesamte Durchbiegung der Stabachse in z-Richtung 
w= = — dimensionslose Anfangsdurchbiegung 
Wes i — dimensionslose Durchbiegung 
é — Dehnung in der 2-Richtung im Beton 
a2 - — Dehnungsgeschwindigkeit 
é- — Dehnung der Bewehrung 


t — Zeit 

o — Spannung in der x-Richtung im Beton; Zug positiv 
P — auBere Kraft 
iP 
Le 


C= — fiktive Spannung 
op — Absolutwert der Prismenfestigkeit fiir Beton 
os — FlieBgrenze der Stahlbewehrung 
o, — Knickspannung der Stahlbetonsaule 
o, — Vorspannung der Bewehrung 
E,x — Materialkonstanten des Betons 


s — #H-10 
E, — Elastizitatsmodul der Stahlbewehrung 


eo — wH|E 
T — X Go 
o, — Spannung in der Stahlbewehrung 
5 a 
a = 
(tf 
» — (4) 


2, — Schlankheit des Stabes fiir die z-Richtung 
A(t) — Kriechfunktion 


Differentialgleichungen des Problems 


Der Stab AB (Abb. 1) ist mit einer Druckkraft P belastet, deren Wirkungslinie 
durch die beiden Gelenke A und B geht. Im unbelasteten Zustand weicht die Stabachse 
von der Geraden ab. Ihre Form ist durch eine Sinusreihe gegeben: 

Wo= >) &m sin m x €. (1) 
m=1 
Es sei gleich hinzugefiigt, daB im Ausdruck (1) die Koeffizienten «, fiir m > 1 klein 
gegeniiber «, vorausgesetzt sind. Dann darf man mit hinreichender Genauigkeit 
schreiben 


Wo = dG Sin 7% é. (2) 


Die Gleichgewichtsbedingungen der Krafte in der x-Richtung und der Momente 
um die y-Achse lauten: 


(cae a FH, Tle), (Gee nee (3) 


va 


[zodF+4 uP Bye [(e)ene — (€) 


F 


ee |e (4) 


z 
& 
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Die Eisenbewehrung habe vor der Betonierung eine Vorspannung o,. Dann gilt 
die Beziehung 

Ee = E+ Ee (5) 


‘ a Bernoullische Hypothese des Ebenbleibens der Querschnitte wird ausgedriickt 
ure 


e=(e),-.— (+) ¢ 24S. (6) 
le 49 — 
saa 43+ 
Z 1 Hj 
Ly 
aa 
ZW 
Pia it es a b-200M 
@- 120M 
Abb. 1 


Damit lauten jetzt die Gleichgewichtsbedingungen ( 3) und (4) 


i 
1 v 
x J oat +ub| @):0+ $ |= —oo (7) 
—1 
5 2 
1 -e\% 02 (w— w = 
x | oodt—pE.(7) Ae aw. (8) 
=} 


Zur Beschreibung der Spannungen in Abhangigkeit von Ort und Zeit wird die 
Zulassigkeit der folgenden Reihenentwicklung vorausgesetzt: 


a (€,¢,t)=— oa a Gye Dow (t) sin ma eae} Sie (t) P,(¢)sinma & (9) 
is Ge Be ee hail 
Hierin sind P, (¢) die Legendreschen Polynome", 
P, =] 
P, (¢) == G 
P, (¢)= 3 (36? —1) 
Ps (C)=3 (508 —3¢) 


u Magnus und Oberhettinger: Formeln und Satze fir die speziellen Funktionen der mathe- 
matischen Physik. Berlin: 1948, S. 70. 
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Analog zum Ausdruck (9) werden fiir die gesamte Durchbiegung des belasteten 
Stabes und fiir die Dehnung in der Stabachse folgende Ansatze verwendet: 


w= >) fm (t)sinmaé (10) 
m=1,2.. 
und 
(e),-0 = 9 (t) + D> em(t)sinma &. (11) 


i = Ny Ce 
Die Gleichgewichtsbedingungen bestimmen die Beziehungen zwischen den Koeffi- 
zienten der Ausdriicke (10) und (11) und des Ausdruckes (9). Fiir die Verhaltnisse der 
Festwerte werden dimensionslose Parameter e, t, s und 6 eingefiihrt. Damit ist, wenn 
man_ bezeichnet 


Dis ee 3 
& neg lO (2) 
€m = €m 10 ; 
3 T 
6. == a 1 — 6) | 
0 08 (Qo (13) 
| 
Cn = —_ Om 
os 
und 
1/3 Ami + (g/t) y sme m? dm 
tn = = ee (14) 


Wie in der Einleitung schon festgestellt wurde, hat das Spannungs-Dehnungs- 
Gesetz die Form: 


=f #0) + A(t) 


(15) 
Die Ausdriicke (9)—(14) und die Beziehung (6) werden in das Spannungs- 


Dehnungs-Gesetz eingefiihrt. Daraus erhalt man die folgende Gleichung: 


: Vy: . 24 inmzé 
Bigg - Daas me t\'— 4 pate aes 
es m m 1 — (9/t) y 8 2 m? 


4 alee Z 6h Se : 2 s 
+ * |e + a An SIM 7 E 2; 3 Om n Py (¢) sin m x :| mie *| oto + 
CPE m m nN 2 

eal A 4 / J . . 

ate =D) (A Gm) SInma & + > : (Ay Qmn)° Pn(¢) sin ma € + (16) 
m m n 

4S) SD) 6,acsin ma ésinkare ee D> > > (Qmn@) P,snma ésinka &é+ 

m k m n k 

Dm end, kde Pe Pes ma gsinkaé| + 

m n k 1 


+ A(t)> |ao+ 2" sin ma & + LMS nce sii ma §| = 05 
§ m m n 

Nach dem Galerkinschen Verfahren wird jetzt diese Gleichung mit sin iz & 
multipliziert und tiber € von =0 bis €=1 integriert. Es entsteht ein unendliches 
System von Gleichungen. Die Gleichung fiir ein gewisses m=7 hat die Form: 
4 ae Nl eee) : 1 Lr . 
Fae |( 2 +40) tn +4 (0) a4] + 5 = [ (FP +20) di + 

is ey atid are 


Sint eae T Ap) Gi ate (1 + FO) 2) di, Py (2 +1 ig) &' ain Pa] + 


n 


lel yyy: 
+45(- 
4m 8 


aa es Ln (4% + Qk n P,) YP (m, k, 2) DD aI PE Gas (a4 -+ P; aki) Y (m, k,i)] xe 
Mm tk en EE ens: 


1 px @? ait Px -(0 
2 lovee? 


(17) 
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wo 
Beh agg nth Os _ taste 
4 1— cos(m+k—i)a 
m+k—it 
ist. 


In der weiteren Rechnung werden alle hodheren Harmonischen in der Sinusreihe 
vernachlassigt. Diese Vernachlassigung ist bis zu einem gewissen Grade durch die 
Voraussetzung begriindet, das in dem Ausdruck (1) fiir die Anfangsdurchbiegung der 
Stabachse die Glieder, die den héheren Harmonischen entsprechen, klein sind. Damit 
bekommt man fiir 1=1 eine Gleichung mit oo! unbekannten Funktionen der Zeit: 
GMT Ot iaklyas 305, OT ASS 835 


5 (2 +144) d+ h(t) a : a |(E* + ra) ey + 


8 


(A agua, = (1-Uea,) Das, Poe (A ta.) Dy an P,| = 


ve al a [se Bees es 


EP 0,4 (ay 4 Py Os )|=0. 


Der Ausdruck in der ersten Klammer ist Null, da er das Spannungs-Dehnungs- 
Gesetz an der Stelle €=€—0 darstellt: 


¢ 0 1 
[é 7. (1—x0) A(t) | =e 5 lo-# [#2 Tf + 1 dg) Gy +A (ta = 0. 
g=t-0 
ieee Gleichung wird noch folgenderweise umgeschrieben : 


ao 


SE ts 


i 


A (t). (19) 

Im weiteren Rechnungsgang werden alle Glieder in Gl. (9), die mit Legendreschen 
Polynomen der Ordnung n > 3 multipliziert sind, vernachlassigt. Es wird also voraus- 
gesetzt, dafi diese Glieder die Form der Spannungsverteilung im Querschnitt nur 
sehr wenig beeinflussen. Nach Multiplikation mit Py bzw. P, (¢), Ps (¢), Ps (¢) und — 
Ausniitzung der Orthogonalitatseigenschaften der Legendre-Polynome erhalt man 
ein System von vier nichtlinearen Differentialgleichungen in a, (¢), a,, (¢), a, (¢) und 
3 (t). Dabei treten folgende bestimmte Integrale auf: 


| aa firs OPQdt=— | POPOPdt= x 


pPredt=a° | P.O) PAO dt= aR. 


aT = 
Das System der nichtlinearen Differentialgleichungen wird in einer Form geschrieben, 
welche die Linearitaét dieser Gleichungen in den Ableitungen @,,, 4,2, @,3, @, betont: 


Cy1 Cyn Crp Cy ayy Cy 
Coy Cop Cog Coa Dg eta REE (20) 
C31 C32 C33 Coa Ay C3 
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Die Funktionen c,, bis c,, und c, bis c, sind gegeben durch 


sp 7 1 T tT 472 27? 
Citar ; -— 7 ly — 75 Org — 3 1 
12) 1— (0/t) y + 4 15a Tt 
272 
472 6 7? T 72 ; 8 12  2t 
> —— ‘ -- a 
PU Pea Gear 352 as) Og tg 0 T Bae ee 
6 72 272 
C13 35 7 112 Cru Fa O13 
27? 272 
Ce ae OL CM i UL 
472 202 272 
OE SO a lL Wo 18 Can = 15 x 712 
if 7 472 272 20? 
Cree et ag a, Ot Cas = D1 718 
272 8 72 t lt+o tT 272 
= ens 
C3 = a7 Gut eq ts ae a 4 ON aged 
272 T 
Ca 5 12 C= GLA) +7 4] i 
272 Raines ASLO 
Co by 12 t 
nie Qg= {ls (t) +7 dM] a, 
27? oe. 
(32 5 TT 45x (21) 


In der Form (20) angeschrieben, ist das System der Differentialgleichungen zur 
numerischen Integration mit Hilfe des Runge-Kutta-Verfahrens! geeignet. Dazu 
sind aber noch die Werte der unbekannten Funktionen ad, @;, @12, 3, 4, fir [=0 
erforderlich. 


Die Anfangswertaufgabe 


Bisher wurde noch nicht erdrtert, wie die Belastung auf den Stab aufgebracht 
wird. Es wird nun wie iiblich vorausgesetzt, da der BelastungsprozeB hinreichend 
langsam vor sich geht, damit keine Tragheitskrafte auftreten, jedoch schnell genug, 
daB wahrend der Belastung keine wesentlichen Kriechverformungen stattfinden. Als 
Augenblick {=0 wird jener Zeitpunkt bezeichnet, in dem der BelastungsprozeB eben 
beendet ist. 

Mit der Voraussetzung, dafi wahrend des Belastungsprozesses keine Kriechver- 
formungen auftreten sollen, ist die Beziehung zwischen Spannungen und Dehnungen 
fiir t=0 gegeben durch 


Oo x 
e=g(1— Fo). (22) 


Die Anfangswertaufgabe ist durch die Gleichgewichtsbedingungen (7) und (8), 
die geometrische Bedingung (6) und durch das Spannungs-Dehnungs-Gesetz (22) 
definiert. Mit denselben Ansatzen fiir o, ¢ und w wie friiher, wobei aber bei allen 
Funktionen von ¢ der Wert fiir = 0 zu nehmen ist, bekommt man ein System von vier 
nichtlinearen algebraischen Gleichungen in a, (0), @, (0), @,3 (0), a, (0), das Ahnlich 
gebaut ist wie das Gleichungssystem (20): 


dyy dy. dyg Aya | [41 (0) dy 

doy A dg dog} | Ay2(0) ee ds (23) 
ds de, dss 34] \ @y3(0) ds ]- 2 
dy, 4, 43 Ugg} \A14(0) dy 


© R. Zurmuth1: Praktische Mathematik fiir Ingenieure und Physiker. Berlin: Springer-Verlag, 
1953, S. 339—348. 
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Die Koeffizienten d,, bis d,, und d, bis d, sind dabei durch die nachstehenden Ausdriicke 
gegeben: 


8 p72 1 tT 27? t? 
dy, el 12 1 —s (o/T) Ss y m2 ma 4 — ca Ay (0) ae bx Ay. (0) ci On Ay (0) 
a 272 37? t 72 472 2 
dya= — (Tey tas (0) + ge a,3(0)) dag = q+ 7% (0) + ge Gas (0) + 3% (0) 
3 T 72 
eye Bbq 12 (0) dz4= 35 113 (9) 
2 72 
dy= 35 “1 (9) dy = 9 #11 (0) 
272 TT 2 
do = On M1 (9) + 7x “1s (9) dissents (0) 
es Te 272 72 72 
do. = 4 == 4 % (0) + Da M12 (0) + 35% (DO) 0 == 1x M13 (0) 
ce , Ee ih ae 2 2 
dog = 2 Ay, ( )+ 63 213 (9) doe, aie Fy (0) te (0) 
pt Ge Bb sx cal 
do4 = aq 12 (0) dy = 4 1—(o/t)ysr 
a 
7 472 : 
dg9= Bell (0) + 25 3 (0) oe 


- Der Koeffizient a, (0) ist von den anderen Unbekannten unabhangig und wird als 
Wurzel der quadratischen Gleichung 


++ 1+6 
2+ Fa(0)| =~ (25) 


gewonnen. Gleichung (25) ist der friiher erhaltenen Gl. (19) analog und stellt das Span- 
nungs-Dehnungs-Gesetz fiir den Punkt = ¢€=0 dar. Die gréBere Wurzel dieser Glei- 
chung hat keine praktische Bedeutung. Die kleinere Wurzel ist: 


a (Oe el yiz2cees 9). (26) 


Sie wird in die Gl. (24) eingetragen. Mit diesem Wert ist auch die Funktion ay (t) 
gemaB Gl. (19) bestimmt. 

4, (0), Ayo (0), G13 (0) und a, (0) stellen die Ausgangswerte fiir das Runge-Kutta- 
Verfahren dar. 


ay (0) |= 


Numerische Beispiele 


Die numerische Rechnung wird fiir den Fall einer quadratischen Saule mit der 
Seite 2b =c=40 cm durchgefiihrt. 
Die Materialkonstanten* op, H und x beziehen sich in allen Beispielen auf Beton der- 
selben Qualitat: 
op = 240 kg/cm? 
E = 240 - 108 kg/cm? 
# = 2 (EH ep — op) /op? = 1/120 cm?/kg, 


wobei fiir den Absolutwert der Stauchung ep = 2°/), genommen ist. Damit ist s= 2 
za nehmen. 


* s. FuBnote?®. 
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Die Kriechfunktion 4 (t) wurde nach Kurven, welche die Ergebnisse der Versuche 
von Glanville und Davis darstellen, berechnet*. Diese Ergebnisse besagen, da® die 
Kriechverformungen bei einer konstanten Last linear abhangig von den Spannungen 
sind und mit der Zeit nach einer Exponentialfunktion der Form: 


A (t) = 2,17 e-0542! 


abnehmen. Die Zeit wird hier in Einheiten von 100 Tagen gemessen. 

Die Bewehrungszahl y hat in allen 24 Beispielen den festen Wert «= 0,012. Der 
Abstand des Bewehrungsschwerpunktes von der Querschnittsachse ist mit ¢e= 12 em 
angenommen. Weiters wird os/op=5000/240 und #,/H=10 angenommen. Damit 


wird 9 = 0,12. 
Die Beispiele beziehen sich auf zwei verschiedene Schlankheiten: 
A —— 75 
AE 100 
zwei Belastungsstufen: 
0,4 
ory == eer OK 


Se OLOtor 
und auf zwei Werte fiir die Vorspannung der Stahleinlagen: 


Lp 900 kg/cm? 
= —— 1350 kg/em?. 


Die kritische Spannung o; wird fiir den geraden Stab nach den Tabellen, die fiir 
Stahlbetonsiulen vorliegen, berechnet!*. Dort ist o, in Abhangigkeit von #, und 
“ os/op zu finden. 

Die Durchbiegung in der Mitte der Stabachse fiir den unbelasteten Stab wird 
gleichfalls variert: 


Das letzte Beispiel bezieht sich auf eine unbewehrte Betonsaule. 

Damit sind die Werte der konstanten Parameter, die in die Gleichungssysteme (20) 
und (23) und in die Gl. (19) und (26) eingehen, festgelegt. Die 25 verschiedenen Kombi- 
nationen dieser Parameter sind in den Tab. 1 bzw. 2 gegeben. 


Numerische Auswertung 


Das System der Differentialgleichungen (20) und die Gl. (19) kénnen fiir das 
numerische Verfahren in der Form 


Te ON as meOve <b Nias C4 
OG, Cig. Gg, Oa Ay Cy 
O Ce, Cog Cog Cag Ao J] Co (27) 
O Cs, Csg Cg3 Cg 13 C3 
O. Cy Can Cag Cag ay C4 


* gs. FuRBnote®. 


® K. Jager: Festigkeitsnachweis im Stahlbetonbau. Wien: Manzsche Verlagsbuchhandlung, 
1948, S. 141, 178. 
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Ingenieur-Archiv XTV/2 


Tabelle 1 
Fall N . er) Ov W max 
‘ ie O0/Ok kg/em? | kg/em? | ™ em 
l 0,4 78,32 
2 0,6 117,5 
0,012 , 5 900 Ig 
3 0,4 63,12 ; 
4 0,6 94,68 
5 0,4 78,32 
6 0,6 117,5 
% 0,012 oe ms 500 a 0 | 4 
8 0,6 94,68 Uibtay 
9 0,4 78,32 B 
10 0,6 117,5 
Gl 0,012 Re ns 900 | 03 | 6 is eee 
12 0,6 94,68 a em? 
13 0,4 78,32 ‘ 
14 0,6 117,5 = 5000 2 
15 0,012 0,4 63,12 1350 0,1 2 Os em2 
16 0,6 94,68 
. 0,4 78,32 1 _ 190 | 
0,6 117,5 
20 0,6 94,68 B= 240-10 Bg 
21 0,4 78,32 i 
22 0,6 117,5 
= 0,012 oo) agra (8b 1. 0,8. | 6 
24 0,6 94,68 
‘25 0 0,4 78,32 0 0,1 2 Unbew. Bet. 
Tabelle 2 
Nr. B Y T to) oy 0 
1 0,5333 0,192 0,65266 0,13789 
2 0,5333 0,192 0,97899 0,091927 
3 0,3 0,108 0,526 0,17110 0,10 0,12 
4 0,3 0,108 0,789 0,11407 
5 0,5333 0,192 0,65266 0,13789 
6 0,5333 0,192 0.97899 0,091927 
7 0,3 0,108 0,526 0,17110 0,20 0,12 
8 0,3 0,108 0,789 0,11407 
9 0,5333 0,192 0,65266 0.13789 
10 0,5333 0,192 0,97899 0,091927 
11 0,3 0,108 0,526 0.17110 0,30 0,12 
12 0,3 0,108 0,789 0,11407 
13 0,5333 0,192 0,65266 0,20683 
14 0,5333 0,192 0,97899 0,13789 
15 0,3 0,108 0,526 0,25665 0,10 0,12 
16 0,3 0,108 0,789 0,17110 
17 0,5333 0,192 0,65266 0,20683 
18 0,5333 0,192 0,97899 0,13789 
19 0,3 0,108 0,526 0,25665 0,20 0,12 
20 0,3 0,108 0,789 0,17110 
21 0,5333 0,192 0,65266 0,20683 
92 0,5333 0,192 0,97899 0,13789 
23 0,3 0,108 0,526 0.25665 0,30 0,12 
24 0,3 0,108 0,789 0.17110 
25 0,5333 0,192 0,65266 0,00000 0,10 0,00 
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geschrieben werden. Die Koeffizienten c,,, ..., C44, Cue cr sma. dure die Gl, (21) 
bestimmt und ¢, ist gemaB GI. (19) ; 
ao Z 
C= — A(t). (28) 
120 of @ + Tao 
Die Werte der Funktionen a, (t), a@,; (f), ..., a, (d) fiir t=0 sind aus dem System 


von algebraischen Gl. (23) und Gl. (26) zu ermitteln. 

Wie schon erwahnt, wird die numerische Lisung des Differentialgleichungssystems 
mit Hilfe des Runge-Kutta-Verfahrens durchgefiihrt. Dafiir wird das in Tab. 3 
gegebene Schema beniitzt. 

Das algebraische Gleichungssystem (23) wird mit Hilfe des Iterationsverfahrens 
gelost. Sehr gute Ausgangswerte fiir die Iteration bekommt man durch Streichung 
aller nichtlinearen Glieder in der ersten Gleichung, aller nichtlinearen Glieder auBer 
@,," und ay, a. in der zweiten Gleichung usw. 

Der ganze Rechnungsgang, namlich Iterationsverfahren fiir das System (23), 
Lésung der quadratischen GI. (26) und Runge-Kutta-Verfahren fiir das System (27), 
wurde am III. Institut fiir Mathematik der Technischen Hochschule in Wien program- 
miert und auf der elektronischen Rechenmaschine IBM 650 ausgefiihrt. 


Tabelle -4. Fall Nr. 5 


t [100 Tage] Ag yy Ibe Ay3 ay 
9822408849 6606320349* 4994341048* 6722361147* 405454147* 

1000000050 8296593749 7285551349* | 2551261648* | 9224867046* | 42627667247* 
2000000050 7515139849 7712688749* 1966983048* 3188228045* 4453746847* 
3000000050 7093349249 7972767849* 1791609248* 1854633146 4588987447* 
4000000050 6858515149 8128131949* 1729065948* 1959726646 4676152847* 
5000000050 6725370449 8219910249* 1703616248* 1833788846 4729805347* 
6000000050 6649074649 8273774049* 1692108448* 1709758946 4762026347* 
7000000050 6605083349 8305266349* 1686458148* 1622586546 4781113147* 
8000000050 6579626849 8323638149* 1683511348* 1567189346 4792331947* 
9000000050 6564864949 8334342049* 1681909848* 1533454446 4798896747* 
1000000051 6556294249 8340573749* 1681015948* 1513335646 4802728347* 
1100000051 6551314549 8344200149* 1680508648* 1501468746 4804961347* 
1200000051 6548420149 8346309949 * 1680217848* 1494511446 4806261547* 
1300000051 6546737349 8347537249* 1680050148* 1490446446 4807018247* 
1400000051 6545758849 8348251149* 1679953148* | 1488076046 4807458547* 
1500000051 6545189849 8348666349* 1679896848* 1486695346 4807714647* 
1600000051 6544858949 8348907849* 1679864148* 1485891646 4807863647* 
1700000051 6544666549 8349048249* 1679845148* 1485423946 4807950247* 
1800000051 6544554649 8349129949* 1679834148* 1485151846 4808000647* 


Die Schrittweite h beim Runge-Kutta-Verfahren wurde so gewahlt, daB mindestens 
vier erste Stellen genau sind. Fiir die Falle, wo die Funktionen do, @4,, ..., @;3 und a, 
gegen einen asymptotischen Wert streben, war h = 100 Tage ausreichend. Bei den Fallen, 
wo die Verformungen schnell anwachsen und daher die Lebensdauer kurz ist, muBten 
Schrittweiten h=10 Tage, h=20 Tage, h=25 Tage, h=40 Tage bentitzt werden. 

Die von der Maschine gelieferten Rechnungsergebnisse sind in den Tab. 4 und 6 
fiir die Falle Nr. 5 und Nr. 6 gezeigt. Die daraus ermittelten Spannungen und Durch- 
biegungen sind in den Tab. 5 und 7 angegeben. 

Die Abb. 2 bis 26 zeigen fiir die 25 untersuchten Kombinationen von Belastung, 
Vorspannung und Vorkriimmung (Tab. 1) den nach den Ergebnissen der numerischen 
Rechnungen ermittelten zeitlichen Verlauf der Stahlspannungen o,/cs, der Betonrand- 
spannungen o/o, und der dimensionslosen Durchbiegung in Stabmitte Wmax- Weiters 
enthalten diese Abbildungen die Betonspannungsverteilung in kg/cm? zu verschiedenen 
Zeitpunkten. 
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oe kg/cm? k 2 ke /em?2 om 
Tage tee - gion B/ Wmex 
Druckseite Zugseite ee eee 
0 — 825 540 —125,0 — 20,4 0,2665 
100 — 2587 307 —119,8 — 5,5 0,3410 
200 — 3582 281 —117,4 3,4 0,3882 
300 — 4154 299 —116,2 8,6 0,4169 
600 — 4783 346 —115,2 14,4 0,4499 
900 — 4906 359 —115,0 15,5 0,4565 
1200 — 4930 361 —114,9 15,7 0,4578 
1500 — 4932 366 —114,9 15,8 0,4581 
—a/0 
Tage 
I II TLE. sad] Lier aed V 
0 1,5956 1,3571 1,0032 0,5993 - 0,2610 
100 1,5292 1,2456 0,8381 0,4050 0,0704 
200 1,4987 1,1924 0,7570 0,3018 —0,0439 
300 1,4839 1,1648 0,7117 0,2442 —0,1104 
600 1,4704 151877 0,6686 0,1822 —0,1840 
900 1,4681 1,1327 0,6601 0,1703 —0,1983 
1200 1,4677 1,1317 0,6584 0,1681 —0,2012 
1500 1,4677 1,1316 0,6581 0,1676 —0,2019 
be = 1/3a,, + 2o/ty ma, 
W max Spe ae 19 o/t y me 
(o/o9)1 = — (A F yy + Mz F a3 + 4) 
Vv 
(o/o9)11t = — (A + @, — Mp/2) > 
6/0)1I = — (a + a, F Ay,/V3 + 2V27 ays) 
4 x 38 8B = 
oe = — 10 t (1 +45 — a — a) + ov F Gy 2,4 10 Bu (Wnax — %) cm? ‘ 
Tabelle 6. Fall Nr. 6 
t [100 Tage] Ay Ay ie Ay3 a, 
9296123649 8319928849* 1051257849* 2320956648* 9455636847* 
2000000049 8929616949 8990902549* 1014146449* 1988580648* 1065426148* 
4000000049 8612057749 9634198849* 1008717349* 1804276048* 1190872448* 
6000000049 8335972249 1024651150* 1025780249* 1727974948* 1319649448* 
8000000049 8095201049 1082557450* 1058651449* 1730976348* 1449733748* 
1000000050 7884637249 1137004650* 1102390049* 1792206148* 1579345248* 
Tabelle 7. Fall Nr. 6 
Tage o. kg/cm? kg/em?2 kg/em? ui... 
Druckseite | Zugseite Ta Omax 
—-, | 0 — 2263 701 —196,2 4,7 0,3444 
S| SS 20 — 3062 760 — 199,7 16,2 0,3862 
; a) 40 — 3856 787 — 203,2 27,4 0,4262 
60 —447] 955 — 206,8 38,1 0,4643 
80 — 5090 1077 —210,2 48,2 0,5004 
ap { Tage | — 0/05 ) 
; ‘ I II | III IV V 
Wye 0 1,6702 1,3916 0,9727 0,4488 —0,0400 
i. 20 1,7000 1,3939 0,9331 0,3709 —0,1380 
40 1,7298 1,3986 0,8993 0,3000 — 0,2330 
60 1,7597 1,4052 0,8981 0,2356 — 0,3241 
80 1,7889 1,4133 0,8479 0,1767 —0,4107 
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Ergebnisse der numerischen Beispiele 


Schon auf den ersten Blick ist zu erkennen, da die Beispiele in zwei Gruppen 
geteilt werden kénnen. Bei den ersten ist die Lebensdauer theoretisch unendlich groB. 
Die Durchbiegung wird in diesen Fillen nach gewisser Zeit asymptotisch gegen cinen 
Endwert streben. In den anderen Fallen ist die Lebensdauer begrenzt durch Erreichen 
der FlieBgrenze in den Stahleinlagen und durch Auftreten groBer Zug- und Druck- 
spannungen im Beton. 

Obwohl die Anzahl der numerischen Beispiele fiir sich ziemlich groB ist, ist 
sie doch zu klein, um auf Grund der vorhandenen Ergebnisse die Lebensdauer fiir einen 
beliebigen Belastungsfall und eine beliebige Form der Saule vorausbestimmen zu 
kénnen. Es ist aber méglich, iiber den Einflu8 einiger Parameter auf die Lebensdauer 
der gedriickten Saiule etwas auszusagen. 

Vier Parameter werden in den vorhandenen Beispielen variiert, und zwar die GroBe 
der Anfangsdurchbiegung, «,, das Verhaltnis o9/oz, die Schlankheit und die Intensitat 
der Vorspannung. 

Bei der Belastung o,—0,6 0, ist in allen betrachteten Fallen die Lebensdauer 
begrenzt; sie ist etwas gréBer bei kleiner Anfangsdurchbiegung %,=0,10 und ver- 
schwindet fast bei «, = 0,30. Bei kleinerer Belastung, 0) = 0,4 o;, und kleinerer Anfangs- 
durchbiegung, «,=0,10 bis 0,20, wird in allen Fallen eine Stabilisierung erreicht. Bei 
gréBerer Anfangsdurchbiegung, «= 0,30, wird auch hier die Lebensdauer ziemlich 
verkiirzt. Die Schlankheit der Saule, bei Festhalten anderer Parameter, beeinfluBt 
die GroBe der Durchbiegung betrachtlich. Eine Anderung der Schlankheit von A, = 75 
auf A, = 100 verdoppelt fast die Durchbiegung. Doch die Lebensdauer und die Gréfe 
der Spannungen in der Saule werden dadurch kaum beeinfluBt. Die Vorspannung 
ist mit o, = 900 kg/cm? und o, = 1350 kg/cm? eingefiihrt. In diesen Grenzen ist nur ein 
kleiner KinfluB auf die GréBe der Durchbiegung, auf das Erreichen der FlieBgrenze in 
den Stahleinlagen und das Auftreten von Zugspannungen im Beton festzustellen. 

Bei einer axial gedriickten Saule wird die GroBe der Belastung durch die tiblichen 
Sicherheitsgrade begrenzt. Damit wird die Lebensdauer hauptsachlich von der GroBe 
der Anfangsdurchbiegung abhangig. In dieser Hinsicht sind also gewisse Forderungen 
an die Genauigkeit der Ausfiihrungsarbeiten zu stellen. 


SchluBbemerkungen 


Den gewonnenen Ergebnissen wurde die Giiltigkeit des Hookeschen Gesetzes fiir 
die Stahleinlagen und eine aus Versuchen* entnommene Arbeitslinie fiir Beton zu- 
grunde gelegt. Diese beriicksichtigt das nichtlineare Verhalten des Betons bei Druck- 
belastung sowie das Kriechen, jedoch nicht das lineare Verhalten bei Entlastung und 
das Verhalten in der Zugzone. 

Nun bedeutet allerdings, wie die gerechneten Beispiele zeigen, das Auftreten von 
Zugspannungen im Beton und das Erreichen der FlieBgrenze im Stahl im allgemeinen 
bereits den Beginn des Zusammenbruches der Saule. Nur in manchen Fallen wird die 
FlieBgrenze in den Stahleinlagen schon bei verhiltnismaBig niedrigen Betondruck- 
spannungen erreicht. Um eine Vorstellung von dem Verhalten solcher Saéulen zu ge- 
winnen, ware es daher erforderlich, den Spannungsablauf auch nach Uberschreiten 
der FlieBgrenze in den Stahleinlagen zu verfolgen. 

Weiters tritt in allen Fallen eine Entlastung des Betons, ¢ < 0, in einigen Teilen oder 
sogar im ganzen Querschnitt auf. Fiir eine genauere Rechnung miiBte man deshalb 
das vollstandige Spannungs-Dehnungs-Gesetz fiir den Fall der Entlastung und hin- 


* g, FuBnoten » 19, 
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sichtlich des Auftretens von Zugspannungen beriicksichtigen. Die numerischen Er- 
gebnisse zeigen allerdings, daB die Spannungen in der Entlastungszone verhiltnis- 
maBig klein sind, so da man sich praktisch noch im linearen Teil der Spannungs- 
Dehnungs-Linie befindet. 

Die erhaltenen Ergebnisse liefern zumindest ein qualitatives Bild vom Verhalten 
gedriickter Betonsaulen und zeigen, da& die Gefahr einer kurzen Lebensdauer auch im 
Falle ziemlich niedriger Drucklasten bestehen kann. 


“ M. Zyezkowski: Some problems of creep buckling of homogeneous and nonhomogeneous 
bars. Proceedings of the I. U. T. A. M. Symposium, Non-Homogeneity in Elasticity and Plasticity, 
Edit.: W. Olszak; Pergamon Press, 1958. 

© N. J. Hoff, W. E. Jahsman and W. Nachbar: A study of creep collapse of a long 
circular cylindrical shell under uniform external pressure. Journal of the Aero-Space Sciences, 
Vo]. 26, No. 10, Oct. 1959, p. 663—669. 

(Eingegangen am 18. Médrz 1960) 


Die Stabilitat gekriimmter, liickenloser Eisenbahngleise 
Von E. Engel, Wien 
Mit 10 Textabbildungen 


Zusammenfassung. Mit Hilfe der Energiemethode werden Formeln entwickelt, mit denen die 
Stabilitatsverhéltnisse von Bogengleisen beurteilt werden kénnen. Dabei wird auf das statische 
System des Gleisrostes eigegangen. 


Einleitung 


Die zu erwartenden Temperaturschwankungen erzeugen in durchgehend ver- 
schweiBten Kisenbahnschienen so hohe Krafte, da einerseits Trennbriiche unter 
Zugkraften und andererseits Gleisverwerfungen unter Druckkraften méglich werden. 
Die zahlreichen, bisher bekannt gewordenen Arbeiten iiber die Stabilitatsverhaltnisse 
von Gleisen unter der Einwirkung von Temperaturdruckkraften gehen von der ver- 
einfachenden, jedoch durchaus unzulassigen Annahme aus, daB sich der Gleisrost 
in waagrechter Ebene durch einen elastischen Trager mit konstanter Biegesteifigkeit 
ersetzen laBt. Zahlreiche Biegeversuche von Gleisrosten lassen jedoch erkennen, dal} 
sich die Gleissteifigkeit in waagrechter Ebene in Abhangigkeit von der Belastung, 
der Stiitzweite usw. veraindert, somit ein eindeutiges Tragheitsmoment nicht angegeben 
werden kann. Uber gekriimmte Gleise liegen keine befriedigenden Arbeiten vor. Der 
Verfasser hat bereits in einer Naherungsberechnung!, die sich auf gerade Gleise be- 
schrinkt, das eigentiimliche Verhalten des Gleisrostes bei der Biegung aufgezeigt, das 
nachfolgend genau untersucht wird. 


I. Biegung des Gleisrostes 


_ Durch eine waagrechte Ausbiegung des Gleisrostes erleiden alle seine Konstruktions- 
glieder Verformungen. AuBer den Schienen werden auch die Schwellen gebogen. Ferner 
verdreht sich die Schwellenachse gegeniiber den Schienenachsen, da die Schienen- 
befestigung nachgiebig ist und Schienen- und Schwellenachse nicht in einer Ebene 
liegen. Gleichzeitig kénnen sich die Schwellen gegentiber den Schienen in der Gleis- 
richtung verschieben. Hierin liegt der wesentliche Unterschied gegentiber einem 


1 E. Engel: Lagesicherheit liickenloser Eisenbahngleise. Z-VDI 102, 373 (1960). 
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Vierendeel-Trager, in dessen Ecken der rechte Winkel unveranderlich ist. Das statische 
Modell, das dem Gleisrost entspricht, zeigt Abb. 1. 


YD 


Abb. 1. Schema eines Gleisrostes 
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Der Drehwinkel «, um den sich der rechte Winkel andert, hangt, wie Versuche 
zeigten, linear mit dem Eckmoment M, zusammen 


1 
a= 7M. (1) 


Nach Uberwinden der Haftreibung zwischen dem Schienenfu8 und der Unterlags- 
platte bleibt hingegen das Eckmoment M), = konstant (Abb. 2), sofern der Schienen- 
fu8 nicht eingespannt ist. 

Durch die Eckmomente M, werden die Schienen wellenformig gebogen. Diese 
sekundaren Ausbiegungen iiberlagern sich der primaren Biegelinie y (x). Legt man die 
(bei Vierendeel-Tragern iibliche) vereinfachende Annahme zugrunde, da die Wende- 
punkte der sekundaren Biegelinie jeweils 
in der Mitte zwischen zwei Schwellen 
liegen, dann lat sich der zusatzliche 
Biegewinkel 6 der Schiene an_ einer 


My 
_K 
Ss 
S 
S 
ae) 


BCG hl; 
Drehwinke/ a 


Abb. 2. Verdrehungsdiagramm einer Schiene Abb. 3. Primire und sekundire Ausbiegung 
auf einer Schwelle (schematisch) einer Schiene 


Schwelle (J, = Tragheitsmoment, # = Elastizitatsmodul einer Schiene) nach dem Satz 
von Mohr angeben (Abb. 3) 


é 


P= 3H, M1. (2) 


Auch die Schwellen werden durch die Eckmomente gebogen (Abb. 4). Der Biege- 
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winkel y am Schienenauflager betrigt (Ho = Elastizitétsmodul, Jo = Trigheitsmoment 
einer Schwelle) 


y= (0) syo7, th. (3) 


Durch den Faktor @ (¢) wird der Widerstand der Bettung ausgedriickt, der die Aus- 
biegung behindert (in Abb. 4 strichliert); ¢ ist der elastische Modul der Verschiebung 
einer Schwelle in der Gleisrichtung, bezogen auf die Lingeneinheit der Schwelle. 
In praktischen Fallen wird © (¢) ~ 1, weshalb hier auf eine genaue Berechnung (die 
mit Hilfe der Differentialgleichung der elastischen Linie fiir die Schwelle erfolgen 
kann) verzichtet wird. Die Auflagerkraft Q,,, auf einer Schiene wird 


One ies (8a) 


§ 


wobei ebenfalls y (¢) ~ 1 in praktischen Fallen. 

Zufolge der Langsverschieblichkeit der Schwellen auf den Schienen und im Schotter- 
bett drehen sich die Schwellen. Entsprechend Abb. 5 ruft eine Verschiebung wu, der 
Schwellen im Schotterbett einen Drehwinkel 6 am Schienenauflager hervor 

2Us 


6= Wc). (4) 


8 


Abb. 4. Biegung einer Schwelle. Die strichlierte Abb.5. Verschiebung einer Schwelle im Schotter- 
Linie gilt, wenn man den Bettungswiderstand _bett. Diestrichlierte Linie gilt, wenn die Biegung 
berticksichtigt der Schwelle mitberiicksichtigt wird 


Die Reaktionskraft Q, des Schotterbettes betragt, elastisches Verhalten desselben 
vorausgesetzt, 
2 
(+ 


Die Faktoren VY (¢) und y (€), die neben dem Bettungsmodul ¢ auch die Schwellen- 
steifigkeit enthalten, sind in praktischen Fallen angenahert gleich eins; hier soll auf 
die genaue Ermittlung der Biegelinie einer Schwelle verzichtet werden; sie fiihrt zu 
unitibersichtlichen Ausdriicken. 

Die elastische Schienenbefestigung erméglicht unter der Kinwirkung einer Kraft Q, 
eine Verschiebung der Schwellen auf den Schienen um das Mab wu, 


Qi = "1° UM. (6) 


Die GroBen mu, und My, sowie v, sind Kennwerte der Schienenbefestigung, die durch 
Versuche zu ermitteln sind. 

Die Auflagerkrafte Q, der Schwellen werden auf die Schienen tibertragen und rufen 
in diesen Normalkrafte NV hervor (N > 0 gilt als Zugkraft). Diese Gurtkrafte bewirken 
elastische Dehnungen « der Schienen, so da®B sich die Punkte einer Schiene um das 


Ingenieur-Archiy XIV/2 10 
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MaB uw gegeniiber dem nicht verformten Zustand in der Schienenlangsrichtung ver- 
schieben. Die elastische Bedingung fiir eine Schiene lautet 


du N 7 
Sa, aE (7) 


Die Gleichgewichtsbedingung fiir die Schiene verlangt 


und das Gleichgewicht einer Schwellenhalfte gegen Verschieben erfordert 


Q10=914+ 9:- (9) 


Die geometrische Bedingung in einem Punkt der Schiene lautet 


U + Uy = Us. (10) 
Mit (7) und (8) folgt daraus 
eHF, tt i tee 
Ur hang (11) 
Durch Umformen von (9) mit (3a), (5), (6) und (11) entsteht 
aan Ee ee . =| " 12 
Picea eo (ita) i 


Die GréBe m bedeutet ein stetiges, pro Langeneinheit des Gleises auf die beiden Schienen 
zusammen einwirkendes Drehmoment. Die Beziehung (12) lautet mit den Bezeich- 
nungen 


ees | “| l 
iso) I= GG Ns, age 
in gektirzter Form 
m=—futgu”. (12a) 
Die geometrische Bedingung fiir die Verdrehung einer EKcke 
yi =at+B+y+o (14) 
liefert nach Umformung mit (1), (2), (3) und (4) 
rtd ib 2 F (¢) x 2eP(C)x HF, ,, 
yao mt ae a ! oe (15) 
Dabei bedeutet 
toeen el € s(t) 2H (fd) mie) 
ee (= | LE, | 6HQJa | vy, 8? i (15) 


Die Beziehung (15) driickt aus, daB das riickdrehende Moment m naherungsweise — 
wenn man die elastische Dehnung der Schienen und damit « vernachlassigt — der 
Neigung der Biegelinie proportional ist, vorausgesetzt, daB sich alle Teile des Gleises 
sowie auch die Bettung elastisch verhalten. Mit den Abkiirzungen 


2 ¥(C) x 2me PV (C)xhF 
dat (1 a) =a (17) 
lautet (15) 
] 
y= 7 (m+du+h wu’) (15a) 
und mit (12a) 
! = h aa 
jade Sea as) 


Fir die Ermittlung der Biegelinie y (x) ist die Gleichung der elastischen Linie 
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rile M; 
Ute (19) 


zu verwenden; M; bedeutet das von den Auf eren Lasten im Querschnitt beider 
Sie : é é' 
Schienen verursachte, in waagrechter Ebene wirkende Biegemoment. 


II. Strenge Lésung des Eigenwertproblems fiir das gerade Gleis 


Es wird ein ideal gerades, unbegrenzt langes Gleis im Schotterbett betrachtet, 
das unter der Wirkung einer inneren Langskraft P, bezogen auf die beiden Schienen, 
steht. Sind die beiden Langskrafte der Schienen ungleich groB, dann bedeutet P deren 
Summe. Verschiebt man das Gleis in der Bettung seitlich, so steigt der Bettungs- 
widerstand anfangs (bis zu einigen Zehntel- 
millimetern Verschiebung) linear mit der Ver- 
schiebung an, um sich dann einer oberen 
Grenze p, zu nahern. Das Verformungsgesetz 
ist empirisch zu bestimmen (Abb. 6). Im ela- 
stischen Bereich gilt daher 


p= =e y. 


Im verformten Zustand des Gleises ist das 
Moment in einem Querschnitt, bezogen auf 
beide Schienen, 


we 
oS 
Bertungswiderstand 


seltliche Verschiebung y 


MA = Py" —p—m'. (21) 
hee 3 Nay: : Abb. 6. Der seitliche Bettungswiderstand 
Die Gleichung (19) der elastischen Linie beider (schematisch) 


Schienen lautet sodann 


2S, y™) + Py" +ey—m' =0 (22) 
und mit (12a) und (18) schlieBlich 
(VI 2g 
2BJ,lt*y +|/—2BJ,! =£4(P—p) i= th]w? + 


+|- a eee sin A eee d] uw" — leat (eNO (23) 
Ml lu Me 

Eine Unannehmlichkeit besteht nun darin, die zugehérige charakteristische Gleichung 
sechsten Grades, die lauter gerade Potenzen hat, zu lésen. Dies 148t sich hier umgehen, 
indem man ein Minimalprinzip verwendet, wie es nachfolgend ausgefiihrt ist. Alle 
partikularen Integrale der linearen, homogenen Differentialgleichung (23) haben die 
Gestalt e+‘, Die Randbedingung, daB die Biegelinie in «= + co stetig bleibt, 
laBt sich nur erfiillen, wenn 7,=0. Legt man das Koordinatensystem derart fest, 
da in x=0 auch w=0 wird, dann mu die Lésung von (23) die Form 


u=C sin 7 (24) 
annehmen ; aperiodische Knickfiguren existieren daher nicht. Fithrt man den Ansatz (24) 


ein, so erhalt man wegen C +0 den Eigenwert 


== Min. (25) 


2 
= Gthtj—d 


Die noch unbekannte Wellenlange / ist derart zu bestimmen, dafi der EKigenwert 


10* 
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ein Minimum wird. Dies fiihrt wieder auf eine Gleichung sechsten Grades, jedoch 
ist hier nur eine ihrer Wurzeln von Interesse. 

In praktischen Fallen wird der zweite Ausdruck in der eckigen Klammer von 
(25) oft klein gegen eins. Wird er vernachlassigt, so erhalt man nach Ausfiithrung 


der Extremalrechnung (G=0) die Naherung 


er Ie | 
in 9 fs 
Prept2y2cekHd, 


Es ist zu beachten, daB die Knicklinge /, von der Art der Schienenbefestigung und 
von den Querschwellen, die durch die GréBe «4 zum Ausdruck kommen, in erster Na- 
herung nicht abhangig ist. Die Knicklast ist als Summe zweier Ausdriicke dargestellt; 
daraus ist die Unzulanglichkeit der Ersatztragertheorie fiir den Gleisrost zu erkennen. 
Die Naherung (25a) gilt nicht mehr, wenn die Wellenlainge J, gro8 wird, was jedenfalls 
fiir kleine Werte von c eintritt. Denn ein unendlich langes Gleis, das nicht eingebettet 
ware (C=0, C=x=0), miiBte ,= oo und P,=0 ergeben, was im Widerspruch zu 
(25a) stiinde. In diesem Falle muB die strenge Lésung (25) verwendet werden. 


III. Das Eigenwertproblem des gekriimmten Gleises 


Fir die Untersuchung eignet sich die Gleichung der elastischen Linie nicht, da auf 
das empirische Gesetz des seitlichen Bettungswiderstandes eingegangen werden muB. 
Ks wird daher das Verfahren von Ritz zur Auffindung des kleinsten Eigenwertes 
verwendet. 

Im Hinblick auf (24) und (18) hat die Knicklinie y (2) des geraden Gleises die Form 
einer Sinuslinie. Es ist daher naheliegend, auch fiir das gekriimmte Gleis naherungsweise 
den Ansatz 


y=v+yosin > x (26) 


zugrunde zu legen. Er enthalt die drei unbekannten Parameter v, y, und J. Die Biege- 
momente m, die von den Schwellen auf die Schienen iibertragen werden und die auf 
die Langeneinheit des Gleises und auf beide Schienen bezogen sind, ergeben sich damit 
und auf Grund von (24), (18) und (15a) zu 


12 
m= Varia Yo COS | x. (27) 
pgth)+f—d 


In dem Ansatz (26) bedeutet v eine konstante radiale Verschiebung der Knicklinie 
gegentiber der unverformten Gleichgewichtslage (Abb. 7). Unmittelbar vor dem Knicken 
wird durch die Langskraft P im Gleis aus Gleichgewichtsgriinden bereits ein bestimmter 
Bettungswiderstand p entsprechend p R= P beansprucht. Eine weitere Verschiebung 
nach auBen erfolgt mit einem Tangentenmodul c, der Bettung, wihrend eine Ver- 
schiebung nach innen mit einem Modul der Entlastung c, erfolgt. Es ist stets c, <c, 
und c, ~¢=konst. Wiirde das Knicken von keiner Radialverschiebung nach auBen 
begleitet sein, dann miiBte der Ansatz einer Sinuslinie die Gleichgewichtsbedingung 
gegen seitliches Verschieben einer Welle verletzen, denn die nach auBen und innen 
gerichteten Halbwellen verursachen verschieden gro8e Anderungen der Bettungs- 
widerstande. 
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Die beim Knicken entstehenden Anderungen der potentiellen Energie des Systems, 


das aus dem Gleisrost, dem Schotterbett und starren Widerlagen besteht, sind fiir 
zwei benachbarte Halbwellen: 


1. Die Biegearbeit der beiden Schienen betragt mit (26) 


+1 
AAp=4 (2B J, y"2de= 2 yp. (28) 
=f 


p= 


unvertorim@e Lage 


Ein konstantes Biegemoment in 
den Schienen, das etwa durch 
elastisches Biegen beim Verlegen 
entstehen kann, hat hier keinen 
EHinflu8. 

2. Die Arbeit des auBeren 
Momentes m, das die sekundaren 
Ausbiegungen der Schienen, die 
Verformung der Schwellen und f 2 f ties = 
die Verdrehung der Schienen- Ms 
befestigung und auferdem die Abb. 7. Knickfigur eines gekriimmten Gleises 
Schwellenauflagerkrafte Q, auf 
den Schienen erzeugt, die ihrerseits die Schienen dehnen und zusiatzliche Schienen- 
langskrafte N bewirken, wird mit dem Ausdruck 


+1 41 
A An=—}[m (a+ B+y+)de=—}|my'de 
| =I 
erfaBt. Mit (26) und (27) gibt die Ausrechnung 
Crh 
ea 
TT Se yo? (29) 
per tind 


3. Die innere Arbeit der Langskraft besteht aus einem Anteil, den die Kraft P 
wahrend der Dehnung leistet, und einem anderen, der durch einen Druckabfall ver- 
Y) 
A a hs 
ursacht durch die Radialverschiebung, und aus einem Teil 1 — V l+y?wxty", der 
die Verlangerung gegeniiber der Abszisse ausdriickt, zusammen. Der mittlere Druck- 


abfall ist A P=2 H Fy = wenn man davon absieht, daB der Bettungswiderstand 


gegen Langsverschiebung der Schwellen den Druckausgleich zwischen den Halbwellen 
geringfiigig behindert. Der Energiebetrag ist daher 


ursacht wird. Die Dehnung eines Bogenelementes setzt sich aus einem Teil ver- 


+1 +1 
AAp=+44P | Zdv—P| (2 +4 y")dz= 
=i) =i) 
2 v ade 
=2HF, I> —2P15—5 7 0. (30) 


4. Die Arbeit des seitlichen Bettungswiderstandes wird naherungsweise erfaBt. 
Die Biegelinie (26) wird durch eine andere Annahme ersetzt (Abb. 8). Sie besteht aus 
zwei verschieden groBen, sinusférmigen Halbwellen, die mit stetiger Tangente an- 
einanderschlieBen. 
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Ks gilt demnach 


bi - Le ad 2 if 
Y 
Yi t Yo=2 Yo 
31 
Yr — Y= 21 ee 
iene 
bo ay 
Ya = Yj 50 Pay Im Sinne der Ritzschen Methode erscheint 
Zp Yir Ye sin FE Ly es gerechtfertigt, die angenommene Gestalt 
‘L 


der Biegelinie zu andern, da es lediglich 
Abb. 8. Biegelinie fiir die geniherte Erfassung darauf ankommt, jeweils die Mittelwerte der 
des Bettungswiderstandes einzelnen Energieanteile niherungsweise zu 
erfassen. Die Naherung gemafi Abb. 8 liefert 
fiir v=0 das exakte Ergebnis. Fiir den anderen Grenzfall v= y, hingegen wird die 
Arbeit des Bettungswiderstandes ungenau erfaft; hier wird aber ihr Anteil am ge- 
samten Potential bedeutungslos. 
Die Arbeit des Bettungswiderstandes lautet somit 
+1 lh 0 
A Ap=—p| y dx — te, | ys" da—te,| ys? ee 
SEN] 0 -—h 


=—2p1e—4|a(14 ral f ul =) | ae? (32) 


Die Anderung 4 A der gesamten potentiellen Energie des Systems ist die Differenz 
der inneren und auferen Arbeitsanteile 
AA=AAz +A Ap—(4 An+A A;). (33) 
Sie mu’ in jedem Gleichgewichtszustand einen stationaren Wert hinsichtlich der in 
(33) enthaltenen Parameter annehmen. Die Bedingung, die das Gleichgewicht gegen 
radiale Verschiebung ausdriickt, . 
(A A) LEE Oe 3 EN o\e 
Tima ee ben oe 
dient zur Bestimmung des Parameters v. Fiir groBe Bogenhalbmesser R wird HF, <c, R? 
und man erhalt dann naherungsweise aus (34) 
*¥e, +o, 
Der Ubergang vom stabilen iiber den indifferenten in den labilen Gleichgewichts- 
zustand ist dadurch gekennzeichnet, da das Potential von einem Minimum in ein 
Maximum iibergefiihrt wird. Das Stabilitatskriterium lautet daher 


(34a) 


2 (AA) 
DY" 
und dient zur Ermittlung der Knicklast. Man erhilt 


=0 (35) 


ae d 
= g— 
[22a In BT Ph Uses 
P(l)=pl1—s a Sele 80 
(gt h+f—-d | at (Vo, + Ve,) 
_ 4EF,P [Ve,—Ve, \ 
ap 7 Re bee Ve, = Min. (36) 
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Die Knicklange /, ist wieder derart zu bestimmen, da8 P (J) ein Minimum wird, 
P= min L (hse 

1. Sonderfall. In Bégen mit groBen Halbmessern wird c,=c,=c, und aus (36) 
entsteht die strenge Lésung fiir die Knicklast des geraden Gleises. GroBe Halbmesser 
(in praktischen Fallen etwa R >1000m) haben keinen spiirbaren EinfluB8 auf die 
Knicklast. 

2. Sonderfall. In engen Bégen wird der verfiigbare seitliche Bettungswiderstand p, 
schon aufgebraucht, um das Gleichgewicht aufrechtzuerhalten (p, R= P). Es wird 
dann c,= 0. Dies ergibt aus (36) groBe Knicklangen [vgl. (25a) fiir das gerade Gleis!]. 
Die Extremalbedingung (36) liefert, wenn man den zweiten Ausdruck in der eckigen 
Klammer vernachlassigt, 


1/ Jy ‘ 
wa | ae (36a) 


Da ferner im Hinblick auf (13), (17) und (16) stets f>d ist, naihert sich die Knicklast 
einem Wert 


Prout l2HF, BS, (36) 
der wesentlich kleiner ist als beim geraden Gleis. 


IV. Richtungsfehler und Grenzgleichgewicht 


Ein streng gerades oder kreisformiges Gleis ist praktisch nicht herstellbar und noch 
weniger im Betrieb zu erhalten. Die Gleislage weist Richtungsfehler (Anfangsstérungen) 
auf, die geduldet werden miissen und die an sich das sichere Befahren noch nicht 
gefahrden. Diese anfanglichen Ausbiegungen werden jedoch durch Druckkrafte in 
den Schienen vergréfert. Es liegt somit kein Eigenwertproblem, sondern ein Ver- 
formungsproblem der Theorie II. Ordnung vor. Ware die Gleisbettung und auch der 
Gleisrost unbeschrankt elastisch, dann wiirde sich, solange die Druckkraft P< Py, 
stets eine stabile Gleichgewichtslage einstellen. In Wirklichkeit (vgl. Abb. 6 und 2) 
wird aber mit zunehmender seitlicher Verschiebung der Bettungswiderstand und bei 
bestimmten Oberbauformen auch der Verdrehungswiderstand zwischen Schienen und 
Schwellen schon nach relativ kleinen Verformungen erschépft. Es wird ein Verfor- 
mungszustand erreicht, bei dem das stabile Gleichgewicht in ein labiles tibergeht 
(Grenzgleichgewicht). Die Verformung vergro8ert sich von selbst, es tritt die ,,Gleis- 
verwerfung* ein. Ein stabiles Gleichgewicht wird erst wieder gefunden, bis sich durch 
die Verformung die Druckkraft in den Schienen geniigend entspannt hat. Die Ver- 
haltnisse zeigt der Verlauf der potentiellen Energie des Systems in Abhangigkeit von 
der seitlichen Ausbiegung des Gleises in Abb. 9. 

Folgende Kriterien lassen ein Urteil itber die Lagesicherheit eines Gleises unter 
Temperaturkraften gewinnen: 

a) Die GréBe jener seitlichen Ausbiegung eines urspriinglich geraden oder kreis- 
formigen Gleises, hervorgerufen durch auBere Krafte (Energiezufuhr), die unter einer 
vorhandenen Druckkraft in den Schienen die labile Grenzgleichgewichtslage darstellt. 

-b) Der Einflu8 von Richtungsfehlern, mit denen in praktischen Fallen zu rechnen 
ist, auf die Grenzgleichgewichtslage. 

c) Die GréBe jener Druckkraft in den Schienen, die bei einer vorhandenen Anfangs- 
stérung der Gleislage bereits das Grenzgleichgewicht erzeugt, ohne diese Stérung 
spiirbar zu vergr6Bern. 

Die allgemeine Aufgabestellung lautet demnach: Wie vergr6Bert sich eine vor- 
gegebene Anfangsstérung und welche Gestalt der Anfangsstérung ergibt bei kleinster 
Druckkraft eine labile Grenzgleichgewichtslage mit bestimmter Ausbiegung? Hine 
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genaue Untersuchung wird durch das nichtlineare, empirische Gesetz des Bettungs- 
widerstandes sowie der Verdrehbarkeit der Schienenbefestigung erschwert. Im fol- 
genden wird mit Hilfe der energetischen Methode eine genaherte Untersuchung durch- 
gefiihrt. 

Es wird zugrunde gelegt, daB die Anfangsstérung der Gleisachse im spannungslosen 
Zustand die Form einer Sinuslinie mit der Wellenlinge J und der Amplitude 7, hat 
(Abb. 10). Die Wellenlinge / soll sich bei der weiteren Verformung nicht andern. Der 
Bettungswiderstand p, sei konstant. Diese Idealisierung erscheint angebracht, da 
bereits nach einigen Zehntelmillimetern seitlicher Verschiebung der Bettungswiderstand 
auf seinen Hichstwert angewachsen ist, die in Betracht stehenden Verschiebungen 
jedoch in der GroéBenordnung von Zentimetern liegen. 
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unvertorinte, stabile Lage 


verformrer, stabiler Endzusiand 


Q\ Latspannung ver 
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Abb. 9. Die potentielle Energie bei wachsender Abb. 10. Anfangsstorung, Grenzgleichgewicht 
Ausbiegung des Gleises und Bettungswiderstand eines Bogengleises 


Die Energieanderungen bei der Verformung von der Anfangsst6rung bis zur Grenz- 
gleichgewichtslage fiir zwei benachbarte Halbwellen sind: 


1. Die Biegearbeit der Schienen. Sie betragt analog (28) 


4H 
A Ap= 45") (37) 


2. Die Arbeit des auBeren Momentes m, das von den Schwellen ausgeiibt wird. Sie 
lautet: 
a) Im elastischen Fall (Haftreibung) analog (29) 


at j 
A Ayes al oe 2 

m 91 m2 (Yo "o )s (38) 

a Oe) of aa 

b) im Falle einer Gleistreibung 

2M ‘ M 

i ’ ’ 8 J. 
AAn=— =" | ly’ —n' |\da= —* (Yo—No): (38a) 


—1 
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3. Die innere Arbeit der Normalkraft P beider Schienen. Entsprechend (30) gilt 


2 2P 
A Ap=2EF lay —2 PI — 7 (yp? — 00°). (39) 


4. Die Arbeit des seitlichen Bettungswiderstandes. FaBt man auch hier wieder 
naherungsweise den nach der Bogeninnenseite ausgebogenen Teil der Biegelinie mit 
der Lange /, als Sinuswelle auf, so ist der Arbeitsanteil 


+1 


ie 4 
A Ap=—pp | ly—n | dz= —2pylv——p Ls (Yo— No— »): (40) 
=) 


Die radiale Verschiebung v ist so zu bestimmen, daB alle Krafte in radialer Richtung 
im Gleichgewicht sind. Das erfordert fiir zwei benachbarte Halbwellen 


= =po(1—“). (41) 


Es mu8 daher stets P <p, R. Um die Lange J, zu ermitteln, kann man naherungsweise 
die Biegelinie durch eine quadratische Parabel 


v — (kh) ‘ 
: mom =) oe 
ersetzen. Aus (41) und (42) folgt 


v=|1—(1->3)"| e—m0)- (43) 


Die Anderung des gesamten Potentials erreicht in der labilen Grenzgleichgewichtslage, 
die durch y, gekennzeichnet ist, ein Maximum; die notwendige Bedingung dafiir ist 
a(A A) 
9 Yo 
Daraus gewinnt man die kritische seitliche Ausbiegung y in Abhangigkeit von 
der Druckkraft P, die in den Schienen urspriinglich wirksam war. Das Ergebnis 
lautet mit (37) und (38) 
a) wenn sich die Schienenbefestigung rein elastisch verformt 


to) 
= 9y,(4 de + 4 Ap—A An — A Ay) = 0. 


21P 2 Z 20 HJ,  4HF,1 
; —27F (1 ot) +291 (1— ot + Se!) — no| ew ++ Et + oth] om 
ch ie Qt ES, 4B, 1 
an B ro 1 |P ea+a”aA+ a (a ery] 
b) wenn zwischen Schienen und Schwellen Gleitreibung herrscht, 
2.UP Se 8 M, 20H), 4HF,1 
TP 1 =o) + 2p9b(1—o8 + Zot) + tt yg [PE 4 ST a — ot 
kK 
Yo = IaH, 8 SET nn (45) 
ie be i ae 
Dabei bedeuten die Abkiirzungen 
m2 
P pe 
pal es AS (46) 


a g+h)+f—4 


Insbesondere entsteht die kritische Anfangsstorung, die bereits das Grenzgleichgewicht 
darstellt, durch Spezialisierung von (44) und (45) fiir yy = 79: 
a) Im elastischen Fall 


2p)? 2 
io = oe le— (1—=) ¢°I. (44a) 
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b) Im Falle der Gleitreibung zwischen Schienen und Schwellen ist die Speziali- 
sierung 7)= Y) praktisch nicht verwertbar. Es mu8 in Wirklichkeit stets ein gewisses 
MaB an Verformung vorhanden sein (also y > 7), bis zwischen Schienen und Schwellen 
Gleitreibung einsetzt; vorher gilt das elastische Verhalten, so daB fiir alle praktischen 
Falle (44a) zu verwenden ist! Ansonsten erhielte man zu giinstige Ergebnisse. 

Von Interesse ist der Sonderfall des geraden Gleises (R= co) ohne Antfangs- 
stérung (7)=0). Durch Spezialisierung von (44) und (45) folgt fiir den Fall a) 


yo" = ST: (44b) 


und fiir den Fall b) 


ee ee ee (45b) 


Es sei bemerkt, da wegen der zugrunde liegenden Form der Biegelinie und der 
Idealisierungen namentlich die Energiebetriige des Bettungswiderstandes und im Fallea) 
auch jene der Schienenbefestigung zu gro8 ausfallen, die Formeln (44) und (45) und 
ihre Sonderfalle daher stets zu giinstige Werte liefern. Der Fehler bleibt jedoch annehm- 
bar klein. 

Der Einflu8 der Lange / einer Halbwelle der Anfangsst6rung wird aus dem Spezial- 
fall (44a) ersichtlich. Bei vorgegebener Kraft P wird das kritische Storma8 7,* um so 
grdBer, je langer die Stérungswelle ist. Daraus folgt umgekehrt: Eine vorhandene 
Storung mit der Amplitude 7, erfordert um so geringere kritische Langskrafte P* in 
den Schienen, um das Grenzgleichgewicht zu erzeugen, je kiirzer diese St6rungswelle 
ist. Kurze Storungswellen sind daher gefahrlicher als lange. 

Einen weiteren Anhaltspunkt fiir die Beurteilung der Wellenlange liefert das 
Verhalten der kritischen Druckkraft P* im Grenzgleichgewichtszustand. Man denkt 
sich (44) baw. (45) nach P aufgelést und fragt nach dem Kleinstwert in Abhangigkeit 
von /. Die kleinste kritische Druckkraft P* tritt bei einer bestimmten kritischen Halb- 
wellenlange /* auf, und die beiden kritischen GroBen hangen von der kritischen Ampli- 
tude y,* sowie von der Anfangsstorung 7) ab. Das Ergebnis lat sich nur in dem 
Spezialfall a) fiir R= oo, 7)7=0 geschlossen darstellen. Es entsteht durch Umkehrung 
von (44b) 

Ap? 2B, 
73 yy* ! [2 


P,=p(1+A)4 = Min! 


Dieser Ausdruck erreicht fiir die ungiinstigste Halbwellenlange 


1/2 BJ, yo* 
x 7% 140 
ih =| o 
seinen Kleinstwert 
QEJ,p 
Shee aed) 
Pe uw (1+aA) 4 4y/ nae (47) 
Im Falle b) entsteht analog aus (45b) 
p,— 4Po% | Int BI, | 8 Mal 
6 Yor | 2 Te Yo* 
und 
d P; 0 8 pol 40 Hd; , 8 Mey (48) 


CM RE RUT ip He? 6 yy * 
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zur Bestimmung von J;. Da das Glied, das M,, enthilt, klein ist, kann man in erster 
Naherung 


b=ltAl, Al<l, (49) 
verwenden. Die Korrektur A/ erhalt man aus (49) und (48) 
Alm 2a 1 
e 2 po fs 324 BS, Yo* - (50) 
rt 1? 


a 


Damit lautet die kleinste kritische Druckkraft im Falle b) entsprechend (45b) und 
im Hinblick auf (47) 


8 Mo, 
me Yo* 


/ 2 EJ, po 


Yo 


Pit x 4| ly. (51) 


V. Zahlenbeispiel 


Die Angaben entsprechen dem auf deutschen und 6sterreichischen vollspurigen 
Hauptbahnen verbreitet verwendeten Rippenplattenoberbau K 49; die empirischen 
Konstanten sind grobe, teilweise geschatzte Durchschnittswerte. Die betrachtete 
Schienenbefestigung verhalt sich nach Fall b) des Abschnittes IV. 


s = 150 cm Schienenabstand 
e= 65cm Schwellenabstand 
a = 50 cm Schwellenitiberstand 
Ta LOO 108 kp Dehnsteifigkeit eimer Schiene 
Pg 0-1) (Or kp em: Biegesteifigkeit einer Schiene 
iodo — 3,0 10°kp= cm? Biegesteifigkeit einer (Beton-) Schwelle 
Y= 2° 10? kp/cm Durchschiebemodul einer Schiene 
(aa 6 10 kp: cmijrad Verdrehungsmodul einer Schiene 
€ = 100 kp/cm? Langsverschiebemodul einer Schwelle 
c = 50 kp/cm? Querverschiebemodul des Gleises 
CA ONL) Hk (oC) wo) 
M,,= 3000 kp - cm Drehmoment einer Schiene auf einer Schwelle im 
Gleitbereich 
Py) = 6 kp/cem groBter seitlicher Bettungswiderstand des Gleises. 


Damit erhalt man fiir die Abkiirzungen (5), (16) und (46) 
#= 104 kp/em, w=180000kp und A=0,001. 


1. Die Knicklast und die Knicklinge eines geraden Gleises werden nach (25a) 
1, = 226 cm, 
P,= 180000 + 520000 = 700000 kp fiir beide Schienen. 
_ 2. Die Knicklast und die Knicklange eines Bogengleises mit R= 300m sind nach 
(36a) und (36b) fiir c,=0 
1, = 690 cm, 
P,, = 180000 + 56000 = 236000 kp fiir beide Schienen. 
3. Die gréBte zu erwartende Temperaturerhéhung von +-50° C erzeugt eine Druck- 
kraft in beiden Schienen 


P=2aAtE Ff, =2- 1,15. 10-5 - 50- 1,32 - 108 = 152.000 kp. 
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4. Die Druckkraft P=152000kp erzeugt das Grenzgleichgewicht im geraden 
Gleis nach (44a), z. B. mit einer Halbwellenlinge /=400cm bei einer kritischen 


Anfangsstérung von 
Net = 0; 9 Cr) 


und im Bogengleis k= 300m bei 
net = 0,2,0M. 


5. Das Grenzgleichgewicht wird durch die Druckkraft P= 152000 kp bei einer 
Halbwellenlange 1=430 cm ohne Anfangsstérung durch Einwirkung auBerer Krafte 
quer zum Gleis erreicht, wenn gem (45) im geraden Gleis die Ausbiegung 


Yo" = 2,0 cm 
und im Bogengleis &A= 300 m 
Yor == 0,7 Cis 


6. Eine Anfangsstérung von der GréBe y= 0,5 em verandert die Amplitude der 
labilen Grenzgleichgewichtslage von 5. im geraden Gleis zu 


Ut == 1,0 CM 
und im Bogengleis k= 300m zu 
Yo* = 0,13 cm. 


7. Die ungiinstigste Halbwellenlange fiir das Grenzgleichgewicht ist mit 
P=152000 kp gemaB (46), (49) und (50) in einem geraden Gleis 


1,* = 430 em. 


Wiirde sich hingegen die Schienenbefestigung mit dem Modul u = 180000 kp unbegrenzt 
elastisch verdrehen, dann wiirde die Druckkraft P = 152000 kp in den Schienen ent- 
sprechend (47) selbst bei unendlich groBer Ausbiegung y, noch kein labiles Grenzgleich- 
gewicht erzeugen kéonnen! 


VI. Schlu8folgerungen 


Ein Knicken im Sinne der Mechanik ist, wie das Zahlenbeispiel zeigt, bei geraden, 
liickenlosen Gleisen unter den groBten zu erwartenden Temperaturspannungen nicht 
zu befiirchten. Bogengleise mit klemen Halbmessern hingegen haben nur geringe 
Knicksicherheiten. Die in Versuchen und in der Natur beobachteten Gleisverwer- 
fungen lassen sich im allgemeinen nicht als ein Knicken erkliren. Sie konnen vielmehr 
eintreten, wenn Richtungsfehler oder auBere Krafte das Gleis so weit verformen, 
daB ein labiles Grenzgleichgewicht erreicht wird. Schon verhaltnismaBig geringe 
Storungen, die mannigfache Ursachen haben konnen und die nicht zu vermeiden sind, 
k6nnen die Lagesicherheit bis auf ein geringes MaB herabsetzen. Bogengleise mit sehr 


klemen Halbmessern erscheinen daher — soweit die verfiigbaren Versuchsergebnisse 
eine verlaBliche Auswertung gestatten — fiir eine durchgehende Verschweifung 
fragwiirdig. 


Die entwickelten Formeln lassen die zahlreichen, in Verwendung stehenden Oberbau- 
formen miteinander vergleichen. Die Lagesicherheit wird wesentlich durch die Ver- 
bindung der Schienen mit den Schwellen beeinfluBt. Befestigungsarten, die zwischen 
Schienenfu8 und Schwelle kein Gleiten, sondern nur eine elastische Verformung 
zulassen [Falla)], sind jenen, die bei einem bestimmten Eckmoment ein wenn auch 
begrenztes Gleiten [Fall b)] erméglichen, weit iiberlegen. Eine abschlieBende zahlen- 
mafige Auswertung kann erst durchgefiihrt werden, wenn eine systematische Ermitt- 
lung der bendtigten Kennwerte erfolgt ist. (Eingegangen am 31. Mérz 1960) 
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Eine Analogie in der Membrantheorie der randbelasteten Schale 
Von J. Barta, Budapest 
Mit 1 Textabbildung 


Zusammenfassung. Die vorliegende Arbeit bezieht sich auf Schalen ohne Biegungswiderstand. 
Es wird gezeigt, da die normierten Pucherschen Spannungsfunktionen der randbelasteten 


Schale proportional den Verschiebungskomponenten einer dehnungsfreien kleinen Deformation 
der Mittelfliche sind. 


Summary. This paper is concerned with the membrane forces of thin shells. It will be shown 
that, if the shell is tangentially loaded on its boundary, the normalized stress functions of Pucher 


are proportional to the displacement components of an inextensional small deformation of the 
middle surface. 


Résumé. Le voile mince est chargé & son contour. Les fonctions de tension de Pucher sont 
ici specialisées par une norme. I] sera prouvé que ces fonctions sont proportionelles aux compo- 
santes de déplacement d’une petite déformation inextensible de la surface médiane. 


I. Grundsiatzliche Beziehungen 


x, y, z sind rechtwinklige kartesische Koordinaten im Raume. G (a, y, z)=0 ist 
die Gleichung der Mittelfliche, und sie wird auch in den Formen x = 2 (y, z), y= y (2, 2), 
z2=2z (x, y) angeschrieben. Die Schnittkraftprojektionen 


Qy (Y; 2); Q.(y, 2),  Yz (y, 2); | 
F(z, x), Wpe(@, 0). ee Liesl enn) sect (1) 
Ss (x, y); Sy (x, Yy); Sxy (x, y), | 


die-den Membranspannungszustand charakterisieren, seien in tiblicher Weise definiert. 
Die Werte Q,, Q., Q,. baw. R. R,, Rex baw. S,, S,, Sy sind also die Projektionen 
der Schnittkrafte an einem Schalenelement, dessen Projektion das Element dy dz bzw. 
dzdx baw. dx dy ist, und werden auf die Seitenlingen der Elemente bezogen. Die 
neun Werte (1) in einem Punkte der Mittelflache sind nicht unabhangig voneinander, 
sondern durch die sechs Gleichungen 


R, = (Rex + Sq) 2, | 
0 
S, = (hz al Sy) ae > 
0 
Se (Sey +9) 2) = 9 
i (2) 
Q; = (S.5 +) z) “Oe? 
Oz 
Q: a (Q)z ste Bex) Oy ’ 
d igeyes (Q2+ dies) pe Abb. 1. Die Lage eines Elementes der 


Mittelflache 


miteinander verkniipft, das heiBt, die Gleichungen (2) miissen in jedem Punkte der 
Mittelflache bestehen. Diese Behauptung laBt sich wie folgt bestatigen. Man betrachtet 
das Element PMN der Mittelflache (Abb. 1), dessen Eckpunkte die Koordinaten 


x, Y; zs 

dz 
x+tdzx, Y, z+ —— dx, 
etdx, yt SY dee, Zz 
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haben. Die Summe der nach Richtung x genommenen Kraftkomponenten soll ver- 
schwinden. Diese Forderung kann sowohl in der Form 


— Ry, da —Syy da +R, + dx=0 
als auch in der Form 
—R,,dx—S,, da +8, ay dg=0 


angeschrieben werden. (Die Lastkrafte und die Tragheitskrafte fallen hier als kleine 
GréBen von héherer Ordnung weg.) Damit haben wir die erste und die zweite der 
Gleichungen (2) bestitigt. In ahnlicher Weise kann man die iibrigen der Gleichungen (2) 
bestatigen. 

Die Gleichungen (2) folgen also aus den Bewegungsgleichungen und gelten sowohl 
in der Statik als auch in der Kinetik. Wiirde man den Schnittkraftkomponenten (1) 
beliebige Werte beilegen, so wiirden Q,, Q., Q,. bzw. R,, R., Rz, baw. S,, S,, Sy) im 
allgemeinen verschiedene Membranspannungszustainde der Schale charakterisieren. 
Wenn aber die Schnittkraftkomponenten (1) solche sind, da sie alle einen selben 
Membranspannungszustand charakterisieren, so sind die Bedingungen (2) eo ipso 
erfiillt. 


II. Die Pucherschen Spannungsfunktionen der randbelasteten Schale 


Von nun an beschranken wir uns auf den Fall, wo die Schale sich in Ruhe befindet 
und die Lastkrafte Tangentialkrafte sind, die am Rande der Schale wirken. In diesem 
Falle existieren, wie in der Literatur gezeigt wird, Funktionen 


U (y, 2), V(z,«), W (x, y); (3) 
die den Gleichungen 
ae aU 9 2x US, worl 0 
Oy? dz? OY O02 OY-0g "OR Oye 
ay aV ay a2V ay aV 
Re 9 Ee eNe 
22 eax 9 a2 0x Terex | Ox? de) 0; (4) 
O82 ct Wo ee ees Nd Wien nee aus ly, 0 
Ox? dy? “axdy dxdy | ay? axt ~~? 
aU aU eu 
Gy = Gay Wor tp) (Ue eae 
eV eV OV 
i, — Pant 2 R, =i eae R., == Sean) (5) 
. @&W _ #&w ise WE 
Ss Oye © y Ou wy Ss» ax ay 
geniigen, wobei Q,, Q., ..., Sz, die Bedingungen (2) erfiillen. Umgekehrt charakteri- 
sieren die nach (4) und (5) berechneten Ausdriicke fiir Q,, Q., . . ., Sy), falls sie die Bedin- 


gungen (2) erfiillen, einen Membranspannungszustand, der ein Gleichgewichtszustand 
ist. Die Funktionen (3) sind die Pucherschen Spannungsfunktionen der randbelasteten 
Schale. 

Wir wollen die jetzt erklarten Funktionen (3) wie folgt normieren: Ein Punkt 
wird auf der Mittelflache G (x, y, z) = 0 festgelegt, und es wird verlangt, daB in diesem 
Punkte 

aV aw aw au aU ov 
de ; Oy ox | oe ay ne Ox 


= 0 
sein soll. 
Man kann diese Normierung immer leicht durchfiihren, indem man iiber die 
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linearen Glieder der Funktionen (3) in geeigneter Weise verfiigt, denn die Werte (5) 
werden durch die linearen Glieder der Funktionen (3) gar nicht beeinfluBt. 


III. Geschichtliches 


Die Herleitung der Pucherschen Formeln (4), (5) wurde in den Originalabhand- 
lungen? und in zahlreichen Lehrbiichern? gegeben. 

Kine Spannungsfunktion der randbelasteten Schale hat zuerst M. Lagally? ein- 
gefiihrt. Er bezeichnet sie als ,,charakteristische Funktion‘‘, und die Differential- 
gleichung, aus der die charakteristische Funktion zu ermitteln ist, als ,,charakteristi- 
sche Gleichung. Die Formeln von M. Lagally wurden von W. Giinther* weiter- 
entwickelt. Beide Autoren beniitzten Flachenkoordinaten. M. Lagally® fand eine 
differentialgeometrische Analogie, die besagt, da zwischen seiner Spannungsfunktion 
und einer gewissen differentialgeometrischen Gré8e eine Proportionalitat besteht. 
Andere Analogien tiber die Schalen haben A. L. Goldenweiser®, C. Stanescu und 
V. Visarion® sowie H. Schaefer’ gefunden. 


IV. Der Zweck des Aufsatzes 


In diesem Aufsatze wollen wir zeigen, daB eine Analogie auch fiir die normierten 
Pucherschen Spannungsfunktionen U (y,z), V (z, x), W (x, y) der randbelasteten 
Schale besteht und wie folgt formuliert werden kann: ,,I[st eine kleine Deformation 
der Mittelflache dehnungsfrei und charakterisiert durch die Verschiebungskompo- 
nenten wu, v, w, so gibt es eine tangentielle Randbelastung, die einen Spannungszustand 
hervorruft, dessen Spannungsfunktionen U, V, W proportional den Verschiebungs- 
komponenten w, v, w sind. Umgekehrt erfahrt die Mittelfliche eine dehnungsfreie 
Deformation zufolge kleiner Verschiebungen w, v, w, die proportional den Spannungs- 
funktionen U, V, W der randbelasteten Schale sind.‘‘ (Die Deformation und die Be- 
lastung, die hier erwahnt sind, stehen in keinem kausalen Zusammenhang.) 

Bevor diese Analogie bewiesen wird, wollen wir die Gleichungen (4) umformen. 


V. Umformung der Gleichungen 
Wir behaupten, da die Gleichungen (4) durch die Gleichungen 


E | es) 
0g = “ay |G 3 
ow aU 

i = 6 
|e z |, 0, (6) 
{aU aV 
lor + a2 |o=° 


1 A. Pucher: Dissertation, Technische Hochschule Graz (1931). 
A. Pucher: Uber den Spannungszustand in doppelt gekrimmten Flachen. Beton und 
Eisen 338, 298 (1934). 
A. Pucher: Uber die Spannungsfunktion beliebig gekriimmter diinner Schalen. Proceedings 
of 5th International Congress for Applied Mechanics, Cambridge/Mass. 134 (1939). 
. B. Laffaille: Mémoires sur l’étude générale des surfaces gauches minces. Abhandl. d. Intern. 
Ver. f. Briickenbau und Hochbau, Bd. 3, 295 (1935). 
2 Siehe z. B. das Buch: A. Pfliger: Elementare Schalenstatik, 2. Aufl., 83 (1957). 
-3 M.Lagally: Uber Spannung und elastische Deformation von unebenen Membranen. 
ZAMM 4, 377 (1924). 
4 W. Ginther: Uber das Gleichgewicht an einer randbelasteten Schale. Abbandlungen der 
Braunschweigischen Wissenschaftlichen Gesellschaft 8, 111 (1956). 
5 A.L. Goldenweiser: Teorija uprugih tonkih obolocek. Gostehizdat, Moskva 1953. 
6 ©. Stanescu i V. Visarion: Statiko-geometrieskaja analogija dlja tonkih uprugih 
obolocek s ortotropijej materiala. Revue de Mécanique Appliquée 3, 109 (Bucuresti 1958). 
7 In Kiirze erscheint: H. Schaefer: Die Analogie zwischen den Verschiebungen und den 
Spannungsfunktionen in der Biegetheorie der Zylinderschale. Ingenieur-Archiv. 
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ersetzt werden diirfen, wobei [...]¢ der Wert des Ausdruckes [...] auf der Mittel- 
fliche G (a, y,z)=0 ist. Um diese Behauptung zu beweisen, gentigt es zu zeigen: 

«) daB die im Abschn. II erklarten Funktionen U (y,z), V (z, x), W (a. y) den 
Gleichungen (6) geniigen, 

B) daB jene Funktionen U (y, 2), 
gentigen, auch den Gleichungen (4) se 

y) daB jene Funktionen U (x, y), V (z, z), W (x, y), welche den Gleichungen (6) 
geniigen, vermittels der Formeln (5) solche Werte Q,, Q:, .-., Sx, liefern, welche den 
Gleichungen (2) geniigen. 

Zu «). Werden die Ausdriicke (5) in die erste der Gleichungen (2) eingesetzt, so 
entsteht die Gleichung 


V (z, x), W (a, y), welche den Gleichungen (6) 


aV OA Vi OW \ ox (7) 
Gee lz ou | auay : 


Der aus den im Abschn. II erklarten Funktionen U (y,z), V (z, x), W (x, y) gebildete 
Ausdruck 


aV aw 
fe a ol G 
seiin der Form 
(= (2, &) as OW (x, p) 
dz oy L= 2X (Y; 2) 


angeschrieben. In dieser Form hingt er von den Veranderlichen y und z ab. Seine 
partielle Ableitung nach z la8t sich wie folgt anschreiben: 


eV aV ox PW dx 
Oz" Ozdxn dz | Ox OY dz ~ 


Ks ist also, mit Riicksicht auf (7), 


ae ad ge eh 


Wworaus 


OV (2, x) OW (x, y) a 
( Oz on oy i =f (y) 


folgt. Diese Gleichung wird, gem&8 unserer Schreibweise, in der Form 


aV ow 

let orle =f (8) 
geschrieben. Ein ahnlicher Gedankengang fiihrt zur Gleichung - 

aV ow 

ae + rl = (2). (9) 
Aus (8) und (9) folgt 

av, ow 

| dz | dy le Aees (10) 


wobei @ eine Konstante ist. In ahnlicher Weise gelangt man zu den Gleichungen 


oe tar |. =) (11) 
Ee oz | 
Oy a (12) 


wobei 6 und ¢ Konstanten sind. Infolge der Normierung (siehe den Abschn. II) ist 
a=b=c=0. Damit nehmen die Gleichungen (10, (11), (12) die Form (6) an. Die im 
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Abschn. IT erklarten normierten Pucherschen Spannungsfunktionen U (y,z), V (z, x), 
W (x, y) geniigen also dem Gleichungssystem (6). 
Zu f). Die Gleichungen (6) seien in der Form 


OV ow 

ke 2 ee cB) 
aw ou 

(ee aE Oz iva =e ele) 
0U oV 

( oy dx La =a (ES) 


angeschrieben. Die Differentiation 02/dy dz sei auf (13), d/dz auf (14), 02/dz dx auf 
(14), d/dy auf (15), 0?/dx dy auf (15) angewendet. Auf diese Weise entstehen fiinf 
Gleichungen. Die erste dieser fiinf Gleichungen sei mit —1, die zweite mit 0? a/dy?, 
die dritte mit (@xz/dy)?, die vierte mit 0? x/dz?, die fiinfte mit (@x/dz)? multipliziert, 
und dann sollen die fiinf Gleichungen addiert werden. So entsteht die erste der Glei- 
chungen (4). In ahnlicher Weise ergeben sich die zweite und die dritte. 
Zu y). Wird (13) partiell nach z differenziert, so entsteht die Gleichung 
@2V  8V ox  #&W ox 
oz" dz0n Oz : Ox Oy dz 


=0, 


woraus unter Beriicksichtigung der Formeln (5) die erste der Gleichung (2) folgt. In 
ahnlicher Weise folgen die iibrigen. 


VI. Der Beweis der Analogie 


Die im Abschn. IV formulierte Analogie beweisen wir wie folgt. Wir fassen die 
Mittelflache der Schale ins Auge. Zur Bestimmung der Punkte dieser Flache mogen die 
Koordinaten x, y als Parameter dienen. Die Gleichung G (x, y, z)=0 der Mittelflache 
sei in der Form z=z (x, y) geschrieben. Die Verschiebungen, welche die Deformation 
der Mittelflache charakterisieren, seien in der Form 

U=U (k, y); v=v0 (a, y); w=w (x, y) (16) 
angeschrieben. Eine einfache geometrische Uberlegung zeigt, daB die kleine Defor- 
mation der Mittelfliche dann und nur dann dehnungsfrei ist, wenn 

dz du+dy dv+dzdw=0 
ist. Ausfihrlicher lautet diese Gleichung 
ou ou ov av z 
dix(5 da + $e dy) + dy($o dx Ney dy) +(5- 


\ 


oder, was dasselbe ist, 
Ov OW OZ Ow Oz ow 
(sr + 5p ay) 49 ty + (5 ae ban) ae det 
ou Ov Ow oz Ow OZ 
ap o ee ox Oy a oy 5, ) dx dy=0. 


Diese Gleichung kann fiir beliebige Werte dz, dy dann und nur dann bestehen, wenn 


Ov Ow az 

oy Fir oy oy eal 

ow dz Ou 17 
Ox Ox AE Ox al ( ) 
ou ov Ow 02 Ow dz 


am 2 az 0% oy ary ie 
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ist. Und nun sei statt (16) 
u=uly,2), v=v(z, x), w=w (x, y) (18) 


geschrieben. Infolge dieser Schreibweise nehmen die Gleichungen (17) die Form 


av ow 

Ee 5 eae 

oN ea 19 
ae + a) oe 
ou Ov 

EE sre 7 le = 9 


an. Die kleine Deformation der Mittelflache ist also dann und nur dann dehnungsfrei, 
wenn die Verschiebungskomponenten (18) den Gleichungen (19) geniigen. Diese 
Gleichungen stimmen formal mit den Gleichungen (6) iiberein. Damit ist die Analogie 
bewiesen. (Hingegangen am 5. April 1960) 


Buchbesprechungen 


An Introduction to Advanced Dynamics. Von S. W. McCuskey. Mit zahlreichen Textabbildungen, 

VIII, 263 Seiten. Reading, Massachusetts: Addison-Wesley Publishing Company, Inc. 1959. 

$ 8.50. 

Das Buch ist als Einfihrung in das Gebiet der analytischen Mechanik fir den kinftigen 
Physikstudenten gedacht und in die folgenden Kapitel unterteilt: Grundlagen der Newtonschen 
Dynamik, Hamiltonsches Prinzip und Lagrangesche Gleichungen, Zentralbewegung, Dynamik 
des starren Korpers, Schwingungsbewegungen, Hamiltonsche Gleichungen und Phasenraum, 
Hamilton-Jacobische Gleichung. 

Eine Durchsicht des Buches fordert unwillkirlich zum Vergleich mit dem bekannten und ganz 
ausgezeichneten Lehrbuch ,,Classical Mechanics‘ von S. Goldstein heraus. Vielleicht ist dem 
Verfasser als Ziel eine ,,vereinfachte‘‘ Version dieses Buches vorgeschwebt. Man darf ihm dann 
zugestehen, da er dieses Ziel im groBen und ganzen auch erreicht hat. Eine Ausnahme bilden 
zum Beispiel die wenigen Seiten iiber relativistische Mechanik, welche kaum geeignet sind, das 
Gebiet dem Verstandnis des mit der Materie noch nicht vertrauten Lesers naherzubringen. Dagegen 
sind die Ausfihrungen tiber den Phasenraum (von der Definition der ,,Dichte‘‘ abgesehen) und 
die Wirkungsvariablen gut gelungen. 

Als Vorstufe zum Studium des Goldsteinschen Buches kann das vorliegende Lehrbuch emp- 
fohlen werden. Die Ausstattung und besonders die Figuren sind ausgezeichnet. 

H. Parkus, Wien 


Festschrift Richard Grammel. Zum 70. Geburtstag am 3. Marz 1959. Unter Mitwirkung von 
A. Betz, H. Goriler, K. Klotter, K. Marguerre, B. Mettler herausgegeben von K. v. Sanden. 
Mit 217 Abbildungen, VII, 372 Seiten. (Ingenieur-Archiv: XXVIII. Band.) Berlin-Géttingen- 
Heidelberg: Springer-Verlag. 1959. DM 52.—. 


Im vergangenen Jahr feierte der Altmeister der technischen Mechanik, Prof. Richard Gram- 
mel, seinen 70. Geburtstag. Fiir seine Freunde und Schiler in aller Welt war dies der Anla&, ihre 
Verehrung durch die Uberreichung einer Festschrift zum Ausdruck zu bringen. So vielgestaltig 
wie das Schaffen Grammels ist auch der Inhalt dieser Schrift. Sie enthalt 38 Arbeiten und ist 
damit zu einem vollstaéndigen Jahresband des Ingenieur-Archivs geworden, jener Zeitschrift, 
deren Griinder und Herausgeber Prof. Grammel ist und die heute aus dem wissenschaftlichen 
Leben Deutschlands nicht mehr weggedacht werden kann. Einer der engsten Freunde Grammels, 
Prof. C. B. Biezeno, hat das Geleitwort verfaBt. 

Ks ist natiirlich ganzlich ausgeschlossen, die Beitrige hier im einzelnen zu wirdigen, ebenso 
wie es ungerecht wire, einzelne Arbeiten herauszugreifen. So sei nur gesagt, da& mit dem Heft 
ein ungewohnlich interessanter Querschnitt durch die technische Mechanik entstanden ist, der 
zeigt, welche Fragen derzeit die Forschung bewegen und welch breiten Raum dabei die Behandlung 
nichtlinearer Probleme einnimmt. Jeder, der irgendwie mit Fragen der Ingenieur-Mechanik zu 
tun hat, wird aus dem Heft reiche Anregung schoépfen. H. Parkus, Wien 
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Berechnung von gufeisernen und emaillierten Druckbehiltern. Von G. Matz und P. Gayer. (Ver- 
fahrenstechnik in Einzeldarstellungen: Band 7.) Mit 97 Abbildungen und 11 Aufgaben, XVI, 
233 Seiten. Berlin-Goéttingen-Heidelberg: Springer-Verlag. 1959. Ganzleinen DM 37.50. 


Obwohl das Hauptziel des Buches in der rechnerischen und experimentellen Untersuchung 
des Festigkeitsverhaltens (Beanspruchung, Abplatzen) der Schutzschicht von auf der Innenseite 
emaillierten gufeisernen Druckkesseln besteht, ist es doch fir einen weit gréReren Kreis von 
Fachleuten interessant als nur fiir den Verfahrenstechniker. Denn es beschrinkt sich nicht auf 
eine kritische Gegentiberstellung der speziellen Messungs- und Rechenergebnisse, sondern enthalt 
auch eine ausfithrliche Darstellung der Technik der Me8streifen (strain gages) und eine griindliche 
Sichtung der dem Praktiker derzeit zur Verfiigung stehenden Methoden zur Berechnung von 
Behaltern und Flanschen. Dabei werden auch Warmespannungen mitbericksichtigt und die Frage 
diskutiert, wieweit die Voraussetzungen der Theorie, vor allem also das Hookesche Gesetz, bei 
GuBeisen zutreffen. 

Die Verfasser haben griindliche Arbeit geleistet. Das Buch kann nicht nur jedem im Kessel- 
und Behalterbau Tatigen empfohlen werden, sondern bietet auch dem Elastizitatstheoretiker eine 
der leider so seltenen Moglichkeiten, seine Uberlegungen und Naherungsannahmen im Licht von 
Versuchsergebnissen kritisch zu prifen. H. Parkus, Wien 


Elektronische Reihe, Band 3: HF-Ubertragungsleitungen, 76 Seiten mit 37 Abbildungen, 
DM 6.—; Band 5: Antennen, 79 Seiten mit 55 Abbildungen, DM 6.—; Band 6: Réhren- 
gleichrichter, 78 Seiten mit 31 Abbildungen, DM 6.—. Stuttgart: Verlag Berliner Union. 1959. 


Die ,,Elektronische Reihe‘‘ umfaft eine Anzahl preiswerter Taschenbiicher tiber Grundlagen, 
Bauformen und Anwendungen verschiedener Bauelemente und Schaltungen und stellt eine Uber- 
setzung einer amerikanischen Buchserie dar. Jeder Band ist in sich abgeschlossen und behandelt 
ein kleines Teilgebiet der Nachrichtentechnik, das dem Leser ohne viel theoretischen Aufwand, 
jedoch trotzdem grindlich und klar nahegebracht wird. GroRer Wert wird durchwegs auf eine 
anschauliche Darstellungsweise gelegt, die nicht zuletzt durch Beispiele, Diagramme und Schalt- 
bilder erreicht wird. Jedem Abschnitt sind einige Wiederholungsfragen beigegeben, die das Ver- 
staéndnis des Lesers fir. das Gebotene vertiefen helfen. 

Der Band ,,HF-Ubertragungsleitungen‘‘ befaBt sich, ausgehend von den Grundlagen, mit 
Begriffen, wie Wellenwiderstand, Stehwellenverhaltnis, Leitungsverlusten usw. und schlieBt mit 
einem Kapitel tiber die Anwendungsmoglichkeiten von Leitungen. Hier findet man Wissenswertes 
tuber 4/2- und A/4-Transformationsstiicke, Resonanz- und Stichleitungen. Auch Verzdgerungs- 
leitungen, Leitungsnachbildungen und Hohlleiter werden behandelt, letztere allerdings nur ganz 
am Rande. 

Der Band ,,Antennen“ bringt, auch wieder auf den Grundlagen aufbauend, Erlauterungen 
verschiedener wichtiger Begriffe, wie z. B. Antennengewinn, Richtcharakteristik usw. Ein Kapitel 
aiber Langdrahtantennen und ein anderes tiber Fragen der Speisung und Anpassung von Antennen 
sind wirklich lesenswert. Ein Kapitel tiber die Ausfiihrungsformen von Antennen rundet den Inhalt 
des Bichleins ab, das sicher jedem, der sich mit Antennen befaBt, etwas zu sagen hat. 

Im Band ,,Rohrengleichrichter“ findet man eine Zusammenstellung und Erklarung der ver- 
schiedenen Schaltungsarten fiir Netzgleichrichter. Nicht zuletzt wegen der beigegebenen Tabellen 
und Nomogramme sowie der durchgerechneten Beispiele wird das Buch jedem Praktiker Inter- 
essantes zu bieten haben. 

Jeder Band ist mit einem Stichwortverzeichnis versehen, das die Verwendbarkeit auBerordent- 
lich erhéht. Der Hochschulingenieur wird in der ,,Elektronischen Reihe“ vielleicht allzu sehr 
die verbindende Theorie vermissen und auch sonst niveaumaBig nicht ganz auf seine Rechnung 
kommen. Wer sich jedoch ohne groBe Miihe und ohne mathematischen Aufwand tiber ein bestimmtes 
Sachgebiet informieren will, findet dazu in diesen Banden ein zuverlassiges und geeignetes Hilfs- 
mittel vor. W. Riedler, Wien 


Refined Iterative Methods for Computation of the Solution and the Eigenvalues of Self-Adjoint 
Boundary Value Problems. Von M. Engeli, Th. Ginsburg, H. Rutishauser und HL. Stiefel. Mit 

’ gahlreichen Textabb., 107 8S. (Mitteilungen aus dem Institut fir angewandte Mathematik an 
der ETH Ziirich, Nr. 8). Basel-Stuttgart: Birkhauser Verlag. 1959. Brosch. sfr./DM 17.—. 


Das vorliegende Buch ist das Ergebnis einer fruchtbaren Zusammenarbeit am Institut fir 
angewandte Mathematik an der ETH Ziirich unter Leitung von E. Stiefel. Es zeigt eine wichtige 
und heute leider noch zu wenig ausgeniitzte Anwendungsmoglichkeit elektronischer Rechen- 
anlagen: Die Priifung und der Vergleich von Naherungsverfahren fir gleiche Aufgabengruppen. 
Gegenstand der Untersuchung sind die verschiedenen Iterationsmethoden zur Losung linearer 
Gleichungssysteme im Vergleich zur Gaufschen Eliminationsmethode und die Methoden zur 


Bestimmung von Eigenwerten. 


160 Buchbesprechungen 


Das Buch besteht aus fiinf Abschnitten. Der erste, von E. Stiefel verfaBte, gibt eine Einfiihrung 
in die Problemstellung und setzt die Theorie der selbstadjungierten Randwertprobleme bei par- 
tiellen Differentialgleichungen auseinander, die durch lineare Gleichungssysteme gelost werden 
kénnen. Es wird weiter auseinandergesetzt, wie die Aufbereitung des Problems far die praktische 
Rechnung erfolgt. Im zweiten Abschnitt setzt H. Rutishauser die Theorie der Gradientenmethoden 
auseinander, die durch das Ziricher Institut bedeutsame Férderung erfahren hat. Bei diesen 
Iterationsmethoden werden die Parameter der Iteration von Schritt zu Schritt neu bestimmt, 
so da8 man hier auf Grund taktischer Erwigungen wahrend der Rechnung, die sich formel- 
miGig erfassen lassen, die weitere Vorgangsweise zweckmaBig wahlen kann. Auch eine Strategie 
kann, wie E. Stiefel gezeigt hat, entwickelt werden, so da man von vornherein sozusagen 
einen besten Plan aufstellen kann. Rutishauser hat nun hier den Zusammenhang zwischen 
den einzelnen Gradientenmethoden besonders durchsichtig zusammengestellt. Th. Ginsburg 
berichtet im dritten Abschnitt iiber die Art und Anlage der Experimente, die den Vergleich 
erméglichten. In erster Linie wurde das Plattenproblem, in zweiter Linie das Stabproblem heran- 
gezogen. Im vierten Abschnitt gibt M. Engeli zuniichst die Theorie der ,,Overrelaxation Methods”, 
bei denen die Parameter der Iteration wahrend der ganzen Rechnung konstant bleiben und die 
Wahl des sogenannten Relaxationsfaktors die Giite der Iteration wesentlich beeinfluBt. Sodann 
berichtet er iiber die Ergebnisse der Experimente. 

Am Schlu8 werden die Ergebnisse der Experimente noch einmal iibersichtlich zusammen- 
gestellt und durch entsprechende Tabellen belegte Hinweise fiir die Leistungsfahigkeit der ein- 
zelnen Verfahren gegeben. Es spielt bei der Auswahl des Verfahrens einerseits der bei Verwendung 
einer Rechenmaschine zur Verfiigung stehende Speicherraum und andererseits die Art und Be- 
schaffenheit der Gleichungsmatrix eine Rolle. 

Diese Untersuchungen kénnen in ihrer Bedeutung gar nicht hoch genug geschatzt werden, 
da ihre Ergebnisse es dem, der unter der Fille der zur Verfiigung stehenden Methoden in einem 
konkreten Fall die zweckmafigste Methode auswihlen soll, méglich machen, nicht erst durch eigene 
langwierige Erfahrung sinnvoll zu entscheiden, sondern sich der Erfahrung dieser Arbeitsgruppe 
bedienen zu konnen. E. Bukovies, Wien 


Praktische Spannungsoptik. Von L. Féppl und H. Monch. Zweite, neubearbeitete Auflage. Mit 
163 Abb., 209 8. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1959. Ganzleinen DM 30.—. 


Die zweite Auflage dieses Buches stellt so wie die erste eine Einfiihrung in den Gegenstand 
dar. Den Fortschritten auf diesem Gebiet der experimentellen Spannungsanalyse wurde natirlich, 
soweit dies in den Rahmen dieses Buches fallt, Rechnung getragen. So nehmen z. B. im Kapitel 
iiber ,,Modellwerkstoffe und Modellherstellung‘‘ die in neuerer Zeit verwendeten Werkstoffe eine 
besondere Stellung ein. Auch die Versuchsberichte im dritten Abschnitt dieses Buches sind auf 
den Stand der zur Zeit tblichen Verfahren gebracht, wobei vorwiegend den Fortschritten auf dem 
Gebiet der raumlichen Spannungsoptik Rechnung getragen wurde. AuSerdem wurde dem Buch 
ein neuer Abschnitt beigefiigt, in dem eine Ubersicht tiber die ,,besonderen‘‘ Verfahren der Span- 
nungsoptik gegeben wird. Dies ist sehr zu begriiBen, obwohl die Ausfiihrungen in diesem Abschnitt 
nur referierenden Charakter haben. Immerhin erfaihrt der aufmerksame Leser, da die Ausfiihrungen 
dieses Buches nicht hinreichen, um iiber die spannungsoptischen Verfahren véllig zu informieren. 
Dieser Umstand wird vom Besprecher bedauert, weil es auSer dem Buch von Mesmer (Spannungs- 
optik, Springer Berlin, 1939), welches natirlich im Hinblick auf die experimentelle Praxis schon 
uberholt ist, kee umfassende deutschsprachige Darstellung des Gegenstandes gibt. Die Verfasser 
haben dies aber auch nicht angestrebt, was schon daraus hervorgeht, daB sie beziiglich eines 
systematischen Literaturverzeichnisses auf die Biicher von Coker-Filon und Heywood ver- 
weisen. 

Zur Ableitung der ,,Hauptgleichung der Spannungsoptik’ (Seite 11) ist noch zu bemerken, 
daB (v,—v,)* T, wenn T wie angegeben die Laufzeit der Welle A, ist, den Gangunterschied 
zwischen A, und A, im Augenblick des Austrittes von A, aus dem Modell darstellt. Fiir den span- 
nungsoptischen Effekt ist aber die Gangdifferenz s=c-(T,— 7) der beiden Teilwellen A, und 
A, nach dem Modell mafgebend; dies wurde ja auch von den Autoren zum Ausdruck gebracht 
(vgl. Abb. 9). Bei der Ableitung der Hauptgleichung ware dementsprechend von dieser maSgeb- 


lichen Gangdifferenz auszugehen. Rk. Beer, Wien 
Eigenttimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fiir den Inhalt verant- 
wortlich: Prof. Dr. H. Parkus, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Steyrermitihl, Wien VI, 


GumpendorferstraBe 40—44. 
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Rechnen in der Chemie. Eine Binfihhrung, Von Dr. techn. Ing. Walter Wittenberger, 
4 Offenbach/Main (frither Aussig/Elbe und Bofors/Schweden). 

- 

& 


Im September 1960 erscheint: 
7 _Zweiter Teil: Chemisch-technisches und physikalisch-chemisches Rechnen unter 
Pe Beriicksichtigung der héheren Mathematik. 
. 2 Mit tiber 700 entwickelten Beispielen und Ubungsaufgaben samt Lésungen und 158 Abbildungen. 
4 Zweit er neubearbeitete und ‘wesentlich erweiterte Auflage. X, 498 Seiten. 8°. 1960. 
eo _ Ganzleinen flexibel S 294.—, DM 49.—, sfr. 50.20.—, $ 11.65 


= Das Buch, das auch auf die Anwendung der héheren Mathematik eingeht, soweit diese fiir den Chemiker normaler- 
__Weise in Betracht kommt, wendet sich in gleicher Weise an den Chemiestudenten wie an den Chemiker in der Praxis. 
Durch die Neubearbeitung des vorliegenden 2. Teiles ist erreicht, da8 dieser Band nunmehr einen vdllig selbstaindigen 


a Charakter erhalt. Gegeniiber der ersten Auflage ist die N euauflage in nahezu allen Abschnitten stark erweitert, die 
if b Zahl der Beispiele und Ubungsaufgaben wurde erheblich vermehrt. So wird das Buch nicht nur zum Uben anregen, 
J _ sondern dariiber hinaus das Interesse wecken an der rechnerischen Behandlung physikalisch-chemischer und SES 
-s ‘technischer Probleme. yee 

~~ 


q Friiher erschien: a 
_Erster Teil: Mit 273 entwickelten Ubungsbeispicien: tiber 1400 Ubungsaufgaben samt Lésun- 


gen und 78 oe Vierte, verbesserte Auflage. IX, 330 Seiten. 8°. 1958. 
Ganzleinen flexibel § 108.—, DM 18.—, sfr. 18.40, $ 4.30 


Chemische Laboratoriumstechnik. Ein Hilfsbuch fir Laboranten und Fachschiler. 
Von Dr. techn. Ing. Walter Wittenberger, Offenbach/Main (friher Aussig/Elbe und Bofors/Schwe- 
den). Finfté, verbesserte und vermehrte Auflage. Mit 397 Textabbildungen. XI, 352 Seiten. 
Se et95 7S the ae Ganzleinen flexibel S 132.—, DM 22.—, sfr. 22.50, $ 5.25 


Biss sis AIS Hilfsbuch und Ratgeber vermittelt die erweiterte und modernisierte Set lee wie bisher dem Laboranten und 

‘Fachschiiler durch die Fiille an erprobten Ratschlaigen ein solides, auf die Praxis eingestelltes Wissen. Das Buch ist 

- geeignet, dem Anfiinger tiber die Befangenheit beim Umgang mit Geraten und Apparaturen hinwegzuhelfen und ihn 
mit den Kniffen vertraut zu peches, die oft fiir ein erfolgreiches Experimentieren notwendig sind . 

és ; ‘Chemiker-Zeitung 


Maschinen und Apparate im Chemiebetrieb. Eine Einfihrung. Von Dr. techn. 
Ing. Walter Wittenberger, derzeit Chemiker der AB-Bofors Nobelkrud, Bofors (Schweden), 
 frither Aussig/Elbe. Mit 446 Abbildungen und einer Bestindigkeitstafel der wichtigsten Werk- 
stoffe. XI, 348 Seiten. 8°. 1949. Steif geheftet S 96.—, DM 16.—, sfr. 16.50, $ 3.80 


»»--. Der erste Teil behandelt den Chemiebetrieb im allgemeinen und der zweite die Werkstoffe, und zwar nicht 
nur jene fiir den Bau von Maschinen und Apparaten fiir chemische Zwecke, sondern ziemlich alle, soweit sie fiir tech- 
nische Zwecke in Frage kommen. Der dritte und vierte Teil befassen sich mit dem Maschinen- und dem Apparatebau 
im besonderen. Diese beiden Teile sind sehr ausfiihrlich bearbeitet und vermitteln wertvolle Auskiinfte betreffs Aufbau 
und Behandlung aller Maschinen und Maschinenanlagen, Apparate und tibrige Hinrichtungen fiir Lagern und Foérdern, 
_ Zerkleinern, Sieben und Mischen, Lésen, Kristallisieren und Sublimieren, Verdampfen und Destillieren usw. ... Das 
B Buch ist nicht nur fiir Chemiker aller Richtungen, sondern ebenso fiir Maschineningenieure und Techniker der Maschi- 
- nenbranche sehr zu empfehlen .. .“ Technische Rundschau 


‘ 


‘ Pregl-Roth, Quantitative organische Mikroanalyse. Von Dr. H. Roth, Badi- 
" sche Anilin- und Soda-Fabrik AG., Ludwigshafen a. Rhein, Landwirtschaftliche Versuchs- 
station Limburgerhof/Pfalz. Siebente, vollkommen neu bearbeitete und erweiterte Auflage. 


_ Mit 115 Abbildungen. XIII, 361 Seiten. Gr.-8°. 1958. 
Ganzleinen S 288.—, DM 48.—, sfr. 49.10, $ 11.45 


,,.Das bekannte Standardwerk Pregl-Roth ist in einer neuen, erweiterten und vollstandig tiberarbeiteten Auflage 
erschienen, In der neuen 7. Auflage ist insbesondere das Kapitel ,,Bestimmung von Atomgruppen“ stark erweitert 
und ausgebaut worden. Die bewahrte Anlage des Buches (Allgemeines, Bestimmung der Elemente, Bestimmung von 
Atomgruppen, Bestimmung physikalischer Konstanten) ist beibehalten worden. Das Buch ist fiir jeden organischen 

, -Chemiker, der sich mit der Analyse oder Konstitutionsaufklarung organischer Substanzen zu befassen hat, eine 
: unentbehrliche Hilfe.“ Klinische Wochenschrift 


Zu beziehen durch Ihre Buchhandlung 
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Rahmentragwerke und Durchlauftrager. Von Dr. Ing. habil. ‘Richard “‘Guldanyt 
weiland o. Professor an der Technischen Hochschule Hannover. Sechste, erweiterte Auflage. 
Aus dem Nachla8 des Verfassers herausgegeben und bearbeitet von Dr. Ing. Horst Reim ann, ie °} 
wissenschaftlicher Assistent an der Technischen Hochschule Hannover. Mit 679 Textabbiles 1 ae 
dungen, 96 Tafeln und 34 Zahlenbeispielen. XXIII, 501 Seiten. Gr.-8°. 1959. rae 
Ganzleinen S 504.—, DM-84.—, sfr. 86. 8 $ 20. i 
» Dieses bekannte Werk erscheint nun in der 6. Auflage, die infolge des Todes des Verfassers von Dr. H . Reimann in Me: ' a 
einigen, das Wesentliche der bisherigen Auflagen jedoch keineswegs beriihrenden Kapiteln iiberarbeitet wurde, gna 
nachdem es seit seiner 2. Auflage stets unverindert herausgegeben worden war. Dieses, das _,,Drehwinkelverfahren“ . 
behandelnde Standard-Werk diirfte in allen Fachkreisen so allgemein bekannt sein und in keinem Biicherbrett. eines As 
Statikers mehr fehlen, da8 es einer speziellen Empfehlung sicher nicht mehr bedarf und ein niheres Eingehen auf seinen 


Inhalt tiberfliissig ist.“ pager Schweizerische Bauzeitung— “ i. 3 : 
A is  aaever oP 


Flach evtiadwarka: Einfihrung in die Elastostatik der ater Platten, Schalen tind 
Faltwerke. Von Dipl.-Ing. Dr. h. ¢. Dr. techn. Karl Girkmann, em. ord. Professor an der Tech- 
nischen Hochschule in Wien, wirkl. Mitglied der Osterr. Akademie der Wissenschaften. Finfte, 
verbesserte und vermehrte Auflage. Mit 318 Textabbildungen. XX XI, 632 Seiten. Gr.-8°. 1959. — 

Ganzleinen § 432.—, DM 72.—, sfr. 73. 70, $ 17. 15s 


3 we 
»,Das heute weit verbreitete Standard-Werk iiber Flachentragwerke erscheint in einer fiinften Auflage, welche sich ey 
von der vorangegangenen vierten nur wenig unterscheidet.. . Umfangreiche Erganzungen zum. Text und zu den Litera- 
turangaben erfassen die Entwicklung seit Erscheinen der letzten Auflage.. . Bemerkenswert ist die beispielhafte | Klar 
heit und Verstindlichkeit des Buches, obgleich bei den Ableitungen dak Schwergewicht auf der allgemeingilltigen 
mathematischen Seite liegt. Diese Klarheit und die den Leser ansprechende Verbundenheit mit der Praxis sind wohl 
der Grund fiir die rasche Verbreitung des Werkes. Es bringt wenig fertige Losungen, weist aber dem selbstandig den: 
kenden Ingenieur und Studierenden — dem das Buch bestens empfohlen werden kann — den Weg zur Lésung der Be 
sich stellenden Probleme. _ hale Schweizerische Bauzettung 


Sy 


Die Cross-Methode und ihre praktische Anwendung. Von Dr. Ing. habil. 
Richard Guldan, 0. Professor an der Technischen Hochschule Hannover. Mit 800 ager 
dungen, 75 Tafeln und 46 Zahlenbeispielen. XIX, 472 Seiten. Gr.-8°. 1955. rae tw 

Ganzleinen § 432.—, DM. 72. =. 24 73. 70, s 17. 15 


»,Der Verfasser bringt eine zusammenfassende Darstellung des Cross-Verfahrens... Im ersten Teil sind die stateceds 
Berechnungsgrundlagen fiir Tragwerke mit und ohne Vouten sowie fiir Rahmen mit ‘und ohne Stabdrehwinkel zusammen- 
gefaBt. Der zweite Teil enthalt Berechnungen mit Zahlenbeispielen, die sich an die im Skelettbau tiblichen Ausfiihrungen — 
halten. Die Zahlen- und Kurventafeln im dritten Teil des Buches werden zur Berechnung der Tragwerke mit und ohne © 
Vouten in zweckmaBiger Weise herangezogen... Das Buch ist fiir Praktiker und Studierende des Ingenieurbaues, _ 
die mit dem Cross-Verfahren arbeiten wollen, eine wertvolle Hilfe.*« ¥ VDI Pigpsilelth 


Hochschule in Wien. Mit 34 Textabbildungen. V, 166 saith an 8°, 1959. a 
Ganzleinen § 228.—, DM 38.—, sfr. 38. 90, $ 9. 08 


.- Mit dem vorliegenden Buch ist innerhalb der Theorie der Warmespannungen die Jetzte wosenthebe Liicke ge- a 
sah loasen worden. Die Vielzah] der behandelten Probleme, die durch ausfiihrliche Darlegungen der Grundlagen vor- 
bereitet werden, und ein umfangreiches Literaturverzeichnis erméglichen es dem Leser, sich uber den letzten Stand © 
der theoretischen Erkenntnisse auf diesem Gebiet ins Bild zu setzen. Der theoretisch interessierte Ingenieur wird ; 
hieraus wertvolle ante unsen bei der Bearbeitung noch ungeloster Hinzelprobleme empfangen kénnen. ee 

Der Ratinghtied 


Die Hochspannungs-Freileitungen. Von Dr. Ing. Karl Girkmann, 0. Prof. ant ae y 
Technischen Hochschule in Wien, und Dr. Ing. Erwin Kénigshofer, Obering. der Osterr. Elektri- ” 
zitatswirtschafts A. G. (Verbundgesellschaft) in Wien. Zweite, erweiterte Auflage. Mit 592 Ab- 
bildungen im Text und 124 Zahlentafeln. XV, 655 Seiten. Gr.-8°. 1952. 

Ganzleinen S 564.—, DM 94.—, sfr. 96.—, $ 22.40 
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